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ABSTRACT 
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An example of a minimal complexity simulation helicopter 
math model is presented. Motivating factors are the 
computational delays, cost, and inflexibility of the very 
sophisticated math models now in common use. A helicopter model 
form is given which addresses each of these factors and provides 
better engineering understanding of the specific handling 
qualities features which are apparent to the simulator pilot. 

The technical approach begins with specification of features 
which are to be modeled followed by a build-up of individual 
vehicle components and definition of equations. Model matching 
and estimation procedures are given which enable the modeling of 
specific helicopters from basic data sources such as flight 
manuals. Checkout procedures are given which provide for total 
model validation. A number of possible model extensions and 
refinements are discussed. Math model computer programs are 
defined and listed. 
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MINIMUM-COMPLEXITY HELICOPTER SIMULATION MATH MODEL PROGRAM 


I . Introduction 

A. Background 

The past decade has seen a trend toward increasingly com- 
plex simulator math models. Part of this has been a result of 
more flight control system sophistication and attention toward a 
number of aerodynamic factors, including interactive aerodynamics 
and aeroelastic effects. Another reason is the availablity of 
large, high speed mainframe and mini-computers. Some simulation 
uses such as aircraft design or failure analysis do justify at- 
tention to detail. Other applications, including may handling 
qualities evaluations, may be better served with lesser sophis- 
tication. Since high complexity also carries the burden of high 
cost of engineering labor and computer facilities, one should ex- 
ercise judgment in math model design. Engineering management 
should be concerned when there is a neglect to determine 
precisely the degree of complexity really needed for a given ap- 
plication . 

The purpose of this report is first to discuss the reasons 
for striving for minimal math model complexity and second to of- 
fer an example of a reasonably useful and credible helicopter 
math model form offering real economy in terms of development and 
computational requirements. Evaluation of handling qualities is 
the main application unde: consideration here, but the same kinds 
of factors would apply to other simulation uses. 

The question being considered is really one of math model 
value versus cost. The value must ultimately be expressed as the 
utility of a math model to provide necessary features which can 
be perceived and used by the simulator pilot. One should expect 
that, as a function of complexity, this model utility approaches 
a fairly flat asymptote with some reasonable level of complexity. 
The other side of the coin is the cost of math model development 
and checkout, also a function of complexity. Unfortunately this 
function can be expected to increase exponentially. These con 
trasting relationships are sketched in Figure 1. The obvious 
question for the simulator user is at what level of model com- 
plexity do these two cost/value curves cross. That is, what is 
the point of diminishing returns on model complexity? 



MODEL 

COST 

OR 

VALUE 



MODEL COMPLEXITY 


Figure 1. Tradeoff of Math Model Cost with User Utility. 


Experience typical of that described in References 1 and 
2 has taught that math model complexity alone does not automati- 
cally provide effectiveness in handling qualities simulations. 
Rather, there can be distracting factors which work counter to 
simulation objectives. Ultimately, limited resources prevent one 
from realizing the full potential of an overly complex simulator 
math model. Other limitations can be a lack of flexibility in 
modeling and restricted clarity in the cause and effect relation- 
ships between model parameters and features. These shortcomings 
raise questions about the value of complexity in helicopter math 
models and are a motivation to consider simpler models. 

The following are some of the undesirable effects of ex- 
cessive math model complexity. 

1 . Computational Delays 

Computational lag and delay is a particularly important 
problem resulting from model complexity. As complexity grows, 
computational delay associated with the math model code increases 
and, in turn, compounds overall visual system delay. Computer 
speed is limited by the hardware and software system being used 
and cannot be easily changed. 
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The result of the delays imposed is reduced fidelity. 
NASA Ames, for example, employs both a Xerox Sigma 8 and CDC 
7600 for their Vertical Motion Simulator (VMS). Using the cur- 
rently implemented ARMCOP helicopter math model (Reference 3) 
and the faster of the two computers (CDC 7600), the computationa j 
delay is about 25 milliseconds. The Sigma may require 60 to 75 
milliseconds to cycle. The former speed is acceptable for some 
math model solutions but visual system digital delay of about 100 
milliseconds can still remain. Although the impact of this 
amount of delay has been minimized at Ames (using methods such as 
those presented in Reference 4), more software complexity has 
been added to correct a problem originally caused by complexity. 
It would seem that this is not as cheap and effective as prevent- 
ing the problem by simplifying the model in the first place. 

2. Cost of Resources 

As model complexity and the amount of computer code grows, 
so do the time and effort required to implement, check, and debug 
the code. The time available to do these is often limited and 
can affect overall math model fidelity if neglected. ARMCOP, for 
example, has several thousand lines of code. In checking and 
debugging code in large programs, a certain number of errors will 
go undetected, and the more code there is, the more likely errors 
will persist. 

In addition to checking and debugging code, there is the 
task of determining model parameters needed to represent a 
specific aircraft. The time and effort required to thoroughly 
validate the model against the real aircraft can be an expensive 
part of any simulation. Math models employing look-up tables can 
have hundreds of parameters which need to be set and confirmed. 
Since some of these values are estimates, an iterative process 
may be required. Limits on time and manpower may restrict this 
process and the fidelity of the model. 

Validation of the math model equations (as opposed to math 
model code) is also a process which may require iteration as the 
model is changed. It is possible that errors in the math model 
will exist even as the model is being used in simulation. Again, 
the number of errors which exist and the time required to fix 
them is a function of the complexity of the model. Time and man- 
power restrictions will limit the ability of the users to find 
and correct these errors and thus degrade the fidelity of the 
model. In order to guarantee that a model is completely correct, 
all parts of the model must be exercised. Lookup tables, for ex- 
ample, require that all numbers in the table be verified as well 
as checked for discontinuities. All equations in the model need 
to be checked to ensure they are theoretically sound. With com- 
plex code, it is unlikely that all of the model will be checked 
as thoroughly as necessary and errors can persist in actively 
used models for long periods of time before they are ever noticed 
or corrected. 
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ARMCOP, for example, still exhibits a problem affecting 
maneuvering flight even though the model has seen wide use. This 
involves a large speed loss during sustained turns. Although 
detected, thi3 problem has not been corrected because of insuf- 
f i cent engineering labor resources. Rather problems are "patched 
up" with flight control system modifications (in this case a 
turn-coordinator). Again, complexity is added to fix a problem 
itself arising from model complexity. 

3. Inflexibility 

There is an inherent tradeoff between complexity and 
flexibility in models of dynamic systems. As more components or 
features are added to a model, it becomes increasingly difficult 
and expensive to perform other modifications. One measure of the 
flexibility of a model is its adaptibility to new computer sys- 
tems and languages or to changes in the code. Large sets of code 
are limited to large computer systems. ARMCOP, for example, re- 
quires the use of a mainframe system. In order to work with the 
model, one must have access to such facilities. 

Once code has been implemented on a machine, it must be 
checked and debugged. Modifications for debugging may require 
recompilation. Most such changes are made before the code is 
used for actual simulation, but it is possible that they will be 
made during a simulation. Even simple model changes can consume 
enough time to hamper simulator productivity. It is not uncommon 
for a software modification, followed by a graphical check, to 
require 20 or 30 minutes of simulator occupancy time. 

The ability to add, remove, or modify efficiently the 
dynamic characteristics of a model is another measure of its 
flexibility. It may be desirable, for example, to have a 
helicopter simulation without rotor cross coupling. A model such 
as ARMCOP, in which cross coupling is inherent, does not allow 
easy removal of this feature. In fact, coupling might be 
removed by adding control system functions to suppress the cou- 
pling, thus further increasing the complexity of the model. The 
emphasis in modeling should be on efficiency while maintaining 
adequate fidelity. 

4. Indirectness of Cause and Effect Relationships 

The ability to see the relationship between model 
parameters and model response features is decreased with com- 
plexity. This relationship is important to handling qualities 
simulation work for two reasons. First, is the need to easily 
make changes in model features. Second is the need to trace er- 
rors which appear in the response modes of the model. These are 
fundamental to working effectively with the model. In order to 
modify response features, one must know what parameters are 
responsible for those features and how to change them. In any 
math model, individual parameters tend to become coupled to many 
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features at once making it difficult to change such features in- 
dependently. The more complicated that math model, the more 
impossible is it to manage individual model response features. 

B. Merits of Considering a Simple Math Model Form 

Turning from the above above list of difficulties issuing 
from model complexity, consider some of the direct, positive 
aspects of considering a simple math model form at the outset. 

It would appear that there are compelling benefits for 
general reductions in the levels of complexity exemplified by 
math models such as ARMCOP and GENHEL (Reference 5). This leads 
us to consider ways to find a compromise between math model com- 
plexity and simulator utility. At one extreme are the highly 
complex models which attempt to acheive effectiveness through 
tugh computational fidelity. As mentioned, these models 

encounter practical limits which not only hamper fidelity but 
also reduce their flexibility and clarity between parameters and 
features. At the other extreme are models such as the linearised 
stability derivative form which are easier to manage but which 
may lack fidelity or be restriced to a small operating envelope. 

The merits of a "compromise" model form would thus be cost. 

benefits derived from the achievement of specific 
fidelity features through minimal software program instructions. 

1 . Cost 


, , The cost benefits will accrue through minimizing labor re 
quired to quantify and checkout the math model implementation 
Development of even modest math models typically involve more 
than one man year of labor. If this process can be shortened to 
less than one man-month, the period envisioned for the proposed 
form, then great savings clearly can be realized. 

Simulator math model software checkout can also require 
substantial effort. However, this is often simply limited bv 
time available and the job might not actually be completed prior 
to simulator use. Again the aim is to realize greatlv reduced 
checkout time through software reduction and to make a comprehen 
sive checkout feasible within a short period of time. 

2. Quality 

* The quality benefits come from confidence that specific 

i eatures needed for effective simulation are represented and that. 

COrr S ct u Here Q ualit y arises from the fact that im- 
plementation and checkout tasks which should be done are, in 

3 real sense, quality follows the degree of 
manageability afforded by the simulator software. 

3. Engineering Understanding 
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One of the moat important benefits to be derived from a 
mimimum-complexity math model is in the potential for more 
clearly understanding cause and effect relationships. For ex- 
ample, if a particular hind and amount of cross-coupling is 
desired, then how does one achieve it through adjustment of math 
model parameters? It is possible by having a close, easy-to- 
follow connection between the physical component representation 
and the resulting physical response features. 

Ar important value of engineering understanding is the 
ability to make model adjustments or refinements in a direct, ef- 
ficient a manner as is possible for a physical helicopter model. 

C. Model Attributes to be Considered 

1. Simulator Application 

It should be stressed that in this case the goal of the 
math model is to be an effective tool for simulation. Model 
fidelity alone is not the solution to simulator effectiveness. 
Rather, it is the ability of the model to produce the desired 
results and insights for the given application. Besides having 
adequate fidelity, the model must also be affordable, manageable 
easily modified and checked, and have a reasonably clear cause 
and effect relationship between parameters and response features 
(at least those perceivable by the pilot). 

2. Handling Qualities Application 

Thus we are motivated to turn to a simple model with these 
qualities for helicopter handling qualitttes simulation which can 
be a more effective tool than existing models. Specifically, the 
purpose here is to propose a minimum-complexity model format 
suitable for helicopter handling qualities simulation. 

It should be remembered that many handling qualities in- 
vestigations involve examination of fairly crude and simple 
parameters such as time constants, damping ratios, or static 
gains. Furthermore the precision with which evaluation pilots 
can perceive such changes often can be disappointing to the en- 
gineer. Thus it is not reasonable to expect that high math model 
resolution is really crucial. If a pilot cannot actually observe 
or be influenced by certain math model effects then those effects 
should probably be considered as excessive complication. 
(Unfortunately, there is presently little quantification of just, 
how sensitive a pilot is to various effects, and this is a poten- 
tial application of a minimum-complexity math model.) 

3. Full Flight Envelope Operation 

The model should be nonlinear and apply to the full 
operating range of a real helicopter including rearward as well 
as forward flight, sideward flight, hover, and transition from 
hover to forward flight. Tne model should include at least 


- 6 - 


r 






first-order flapping degrees of freedom and all rigid body de- 
grees of freedom. The higher-order flapping modes and any 
structural modes beyond the frequency range of interest for han- 
dling qualities should not be included unless high-gain flight 
control systems are involvede. 

4. Modularity 

The form of the model will be modular. This will allow 
the flexibility of adding alternative rotor models, if desired, 
as well as other lifting surfaces. Any combination of components 
can be combined including models of pilots and control systems 
making the model adaptable to a variety of helicopters and sub- 
systems. The full utility of the proposed model format will 
become apparant as the structure of the model is described in 
more detail. 

5. Microcomputer Adaptability 

The math model form will be compatible with microcomputer 
use, at least on a non-real-time basis. It has been found that 
math model development and checkout can be done to a large extent 
on small, inexpensive desktop microcomputers. This of course 
demands that the software be reasonably compact. 

D. Report Organization 

The presentation which follows consists of four parts: 
(i) approach to modeling, (ii) matching and estimation proce- 
dures, (iii) checkout procedures, and (iv) extensions and 
modifications of the model. In addition various detailed i ifor- 
mation is contained in appendices. 

1 . Modeling Approach 

In the first section, the modeling approach is described 
in order to establish the theoretical foundation for the model. 
This is useful for understanding, modifying or extending the 
model and for its effective use as a simulator tool. In addi- 

tion, a description of the features and components of this 
specific model is given. The model is used to represent a Bell 
AH-1S Cobra. All parameters and variables from this aircraft are 
provided here along with the actual code. The sample version 
shows the extent of the code in terms of number of parameters, 
number of lines of code, number of computations, etc. and can be 
compared to an ARMCOP version of the same aircraft. 

2. Matching and Estimating Procedures 

In the next section, the matching and estimating proce 
dures used to obtain model parameters are described. The sample 
version of the AH-1S is used as a specific example. The model is 
then exercised and the estimated parameters varied in order to 
tune the model to fit performance data. 



3. Checkout Procedures 


The third section describes several methods of checking 
the math model code. The size of the model and the modular for- 
mat are conducive to efficient checking. Methods are then 
presented for varifying the math model equations and are il- 
lustrated using the sample version. 

4. Model Extensions and Refinements 

Finally, in the last section, possibilities for extending 
or modifying the model are introduced to demonstrate the 
flexibility of the model format. The potential for a much im- 
proved level of simulation effectiveness using these extensions 
and modifications is revealed and explained in terms of the ap- 
proach taken to the modeling process. 
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I I . Technical Approach 

A. Specification of Desired Math Model Features 

The approach to modeling will begin with a list of desired 
features. This list will serve as a specification upon which to 
formulate a minimum-complexity model containing only those com- 
ponents and equations directly responsible for the desired 
features. The model will be customized to the problem being 
studied . 


We shall assume that the model is intended for handling 
qualities simulation and that the features to be included in the 
model should be features which are observable or needed by a 
pilot. It should be assumed that the model will be operated over 
a specified flight envelope and controlled by a given flight con- 
trol configuration. This sets limits on speed, acceleration, and 
frequency response. The features to be included by the following 
example are listed in Table 1. Of course these are subject to 
change depending upon the application. 


Table 1 . Desired Features 


1. First-order flapping dynamics for main rotor (coupled or 
uncoupled ) . 

2. Main rotor induced velocity computation. 

3. All rigid-body degrees of freedom. 

4. Realistic power requirements over desired flight envelope. 

5. Rearward and sideward flight without computational sin- 
gularities . 

6. Hover dynamic modes: 

- -longitudinal and lateral hover cubics 
--rotor-body coupling with flapping 

7. Forward flight dynamic modes: 

--short period and pnuguoid 
--roll mode and Dutch roll 
--rotor-body coupling with flapping 


8 . 

Dihedral effect. 



9. 

Correct transition from 

hover to forward 

flight. 

10 . 

Potential for rotor RPM 

variation . 


11 . 

Correct power-off glide 

for min R/D and 

max glid 


F 
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1. First-Order Flapping 


It has been shown in Reference 6 that rotor flapping can 
couple with rigid-body modes in regions which affect handling 
qualities . This occurs in the lower frequency or ’regressing 
flapping" modes. However, this effect can be modeled with a 
first-order flapping equation in the pitch and roll axes. 

The time constant involved in the regressing flapping mode 
is directly proportional to the product of rotor angular velocity 
and Lock number . Thus only the commonly available rotor mass and 
geometric parameters are needed. 

The actual flapping response is modified by coupling with 
the fuselage at the hub restraint. Since this involves the clas- 
sical rigid body modal reponse, it is discussed further under 
items 6 and 7 below. 

The feature of flapping which is most important to a 
pnot-in-the-lcop simulation is the apparent control lag follow- 
ing cyclic input. This lag is in effect the time required to 
process the tip path plane to a new orientation. A typical value 
for the effective lag is about 0.1 sec — significant because it is 
comparable to the pilot’s own neuromuscular lag. 

2. Main Rotor Induced-Velocity Computation 

A particularly important feature of a helicopter is the 
relationship among thrust, power, and airspeed. This relation- 
ship arises from the induced-velocity of air passing through the 
rotor disc. 

There are a number of complicating factors, but, to a 
first-order approximation, induced-velocity effects can be 
modeled with a classical momentum theory model wherein thrust and 
induced-velocity interact in an aerodynamic feedback loop 
Computation is complicated, however, because this feedback is 
highly nonlinear. 

Another aspect of the induced-velocity is its effect on 
adjacent surfaces. The rotor induced-velocity field impinges on 
the wing, horizontal tail, and fuselage and varies with airspeed 
and flight path direction. 

3. Rigid-Body Degrees of Freedom 

Normally, six rigid-body degrees of freedom are needed for 
useful manned simulation. Pilot workload arises from constant 
attention to roll, pitch, and yaw as well as translation fore- 
and-aft, to the side, and vertically. Only under special 
conditions might one desire to eliminate one of these via. for 
example, the assumption of perfectly coordinated forward flight. 
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4. Power Requirements Over Flight Envelope 

A common source of real aircraft data appropriate for 
verifying a math model is performance data in terms of power re- 
quired for various trim conditions. The power or torque required 
is immediately obvious and important to a pilot and varies sub - 
stantially from hover through transition and finally in forward 
flight. 


Power requirements can be easily computed once main and 
tail rotor induced velocities are established. 

5. Rearward and Sideward Flight 

In a full-flight-envelope model involving circulation 
lifting surfaces, computational singularities can exist, depend- 
ing upon the model form used. These singularities come from 
trigonometric functions for angle of attack, sideslip, etc., but 
are avoided in this model by using a quadratic lift coefficient 
method. For this technique, forces for lifting surfaces are com- 
puted using quadratic coefficients multiplied by the squares of 
velocity components so that negative velocities cannot cause sin- 
gularities. No explicit computation of angle of attack or 
sideslip is needed and, indeed, should be completely avoided. 

6 . Hover Dynamic Modes 

Hovering flight is characterized by similar dynamics in 
each the pitch and roll axes, including sets of high and low fre- 
quency response modes. In addition, the yaw axis contains a 
predominant yaw damping mode. These dynamics can couple with 
regressing flapping dynamics. All are apparent to the pilot in 
operating the aircraft whether trimming, maneuvering, or flying 
unattended . 

Pitch and roll are classically described by the "hover 
cubic," but this generally neglects coupling with the rotor which 
can be important. This is easily computed, however, through in- 
clusion of the flapping dynamics as described earlier. 

The phugoid mode for hover results from the combination of 
dihedral and gravity force. Effective dihedral is particularly 
apparent in unaggressive sideward flight because the pilot must 
continually add lateral control as sideward velocity increases. 

7. Forward-Flight Dynamic Modes 

In forward flight, the dominant rigid body dynamics of a 
helicopter resemble those of a conventional fixed-wing airplane 
and include short-period, phugoid, dutch roll, and spiral modes. 
There is also likely to be significant coupling with flapping 
dynamics . 
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8. Correct Transition from Hover to Forward Flight 

Transition effects are an important part of the piloting 
task when accelerating from hover into the forward flight region. 

These effects are a combined result of a "dihedral effect 
in the x-axis and the varying rotor downwash effect on the 
horizontal tail. 

9. Effects of Rotor RPM Variation 

Rotor RPM can affect helicopter dynamics in a number of 
ways, including thrust, flapping response, and heave damping. 

The effects of rotor speed • • atlon are tied, however, to 

the rotor-engine-governor combination. For a number of applica- 
tions it may be sufficient to assume a constant rotor RPM. This 
will be done here. 

10. Cross Coupling 

A variety of cross coupling effects can be present in 
helicopters. Some of these such as collective-to-yaw coupling 
are easy-to-see, first-order phenomena. These are generally in- 
herent in the basic dynamics if reasonable first-principles 
thrust and rotor models are used. 

Other coupling effects may be more subtle or less predict- 
able and should be added only where needed or desired by the 
simulator user. These can be inserted directly in the equations 
of motion as coupling terms arising from both states and control- 
s. One should distinguish among coupling due to (i) rotor hub 
moments, (ii) flapping dynamics, and (iii) dihedral effects. 

Cross -coupl ing occurs naturally when coupled rotor and 
hub-moment expressions are inciuded. However, these may no suf- 
fice in matching actual cross-coupling observed in a particular 
design. One approach is to begin with decoupled equations then 
systematically add terms which provide a suitable match. (This 
is demonstrated in the A109 example in Appendix D. 

A useful guide to cross-coupling sources is borrowed from 
Reference 7 and shown below in Table 2. 

13. Correct Power-Off Glide 

Helicopters, like fixed-wing aircraft, need to exhibit 
reasonable performance when power is reduced. This can be a 
highly complex issue if ring vortex rotor states are included. 
However, many handling qualities investigations can be conducted 
using only the normal thrust model described above but tailoring 
the full-down collective pitch and aerodynamic drag to yield 
realistic f o -ward-velocity autorotative glide characteristics. 
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Table 2 . 


Sources . 


Single-Rotor Helicopter Coupling 


1 '^^\ s Response 

Input Axis 

kit Ct) 

Ron 

Yaw 

CNmo/Descent 

L ongitodinal stick 

Prime 

( t ) lateral flapping due to 
longitudinal stick 

(2) lateral flapping due to 
pitch rate 

(3) lateral flopping due to 
load factor 

Negligible 

Desired for vertical flight 
path control in forwar J 
flight 

Lateral stick 

( 1) longitudinal flopping 
due to loterol stick 

(2) iongltudmalal flapping 
due tc roll rote 

Prime 

i 

( 1 ) Undsslred in hover, 
caused by directional 
stability 

12) Desired for turn 
coordination and 
heading cont-ol In 
forwardfllgnt 

Descent with bank 
angle at fixed power 

ftuMer 

Negligible 

( 0 Roil due total) rotor 
♦hruat 

(2) Roll due to sideslip 

Prime (hover) 

Undesired due to power 
changes in hover 

Collective 

(1) Transient longitudinal 
flapping with load factor 

(2) Stead/ longitudinal 
flapping due to climb/ 
descent In forward flight 
caused by rotor flapping 

( 3 ) P itch due to changes 
in horizontal tall lift 

( 1) Transient lateral 
flapptngwlth load factor 

( 2) Stood/ lateral 
flapping with climb and 
descent 

(3) Sideslip induced by 
power change causes 
roll due to dihedral 

Povrer change varies 
requirement for tall rotor 
thrust 

Prime 

J 


(Borrowed from Blake and Alansky, AH5 Forum, 1975) 


B. Component Build-Up 


With a specification of desired features, essential model 
components can then be chosen. These components contain the 
mechanisms which provide forces and moments, power dissipation, 
.stability and control, and rotor dynamics. 


The six 
of the above 
•3 lists these 
eaon component 
■his is a list 
starting point 
their phsyical 
below . 


components are considered necessary to provide a l ; 
response features are shown in Figure 2. T n i 
components along with the physical feature:, 
and the resulting response features. F n effect, 

of qualitative model requirements which form a 
for detailed model design. The components and 
elements are described and discussed individually 
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!•) Main rotor: 


2.) Fuselage: 


3 • ) Tail rotor : 


4.) Horizontal tail 


5. ) Wing: 


6 • ' Vertical tail : 


Thrust 

Torque 

Induced velocity 
Tip path plane lag 
Induced power 
Profile power 
No off-axis 

flapping stiffness 
Decoupled TPP dynamics 
Constant RPM 


Mass at C.G. 
Moments of inertia 
Parasite power 
Cross products of 
inertia = 0 


Thrust 

Torque 

Induced velocity 
Induced power 
Profile power 


Lift / Stall 
Exposure to main 
rotor induced vel 


Lift / Stall 
Induced drag 
Induced power 
Exposure to main 
rotor induced vel. 


Lift / Stall 


1st order flap- 
ping 

Power required 
Trim 
Phugoid 
Short period 
Dihedral 
Pitch mode 
Roll mode 
Min x-coupling 
Power off glide 


Trim 

Power required 
Min x-coupling 
Power off glide 


Trim 

Power required 
Roll mode 


Short period 
Trim 

Pitch mode 
Power required 


Trim 

Power required 


Dutch roll 
Roll mode 


1 
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1 . Main rotor 


. The P f imary component of this model is the main rotor. It 

is the main feature responsible for producing characteristics 
unique to a helicopter, in particular, a vertical thrust vector 
and an induced-velocity field. Other key features include rotor 
torque, dihedral effect, flapping stiffness (rate damping), and 
flapping dynamics (tip-path-plane lag). B 

The basis for the model used here is primarily the 
autogyro theory presented by Glauert ir. Reference 8 and extended 

S y e L °5 k v 1 JJu Rer ? re « CO 9 ' The hi3her order flapping dynamics as 
defined by Chen in Reference 6 are simplified according to the 

first order model developed by Curtiss and presented in 
References 1 and 2. 


Thrust and induced velocity are computed assuming a 
uniform flow distribution. As described earlier, the tip-path- 
plane orientation (flapping angles) are modeled as simple first 
order lags giving the main rotor the qualities of a force ac- 
tuator with a lag. The tip-path-plane dynamics can be extended 
using either a coupled first-order model or a coupled second 

order model cased on simplification of Chen’s rotor equations in 
Reference 6 . 


The main rotor model contributes largely to the power re- 
quirement feature of the model. In hover, nearly 80% of total 
power is absorbed by the main rotor, and, in forward flight, it 
is as much as 60%. In hover, rotor downwash on the fuselage also 
contributes to power losses. 


Tip 

rederived in 
References 3 
time hover 
instantaneous 
in body-axis 


path plane and hub moment equations were 
a body-fixed axis system from the equations in 
and 6. This was done in order to avoid the real 
simulation problems which can arise from large 
changes in the wind-axis angles for small changes 
translational velocities. 


f suggested that two major components of cross cou- 
P ing e avoided until the detailed model matching process is 
underway. One of these is the off-axis hub moments due tc flap- 
ping (L al and M bl ), and the second is the off-axis coupling in 


the tip path plane dynamics. It has been found that 
these higher order effects in a simple model does not 
caxly produce a high quality match to flight data. 


including 
automat i - 


dbj /dv 


included through the variables 
the first order flapping equa- 


The dihedral effect is 
and da^/du which appear in 

tions. Values can be can be computed using first-principles 
factors consisting of thrust coefficient and tip velocity The 

J eatur j ,)f responsible for the phugoid-like modes ir 
hover and forward flight. 
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The portion of L bi and M fil due to both hinge offset and 

rotor spring stiffness are included in a separate parameter, 
dL.dAl. Thus, the total flapping stiffness can be directly 
varied through this one parameter. 

Pitch and roll mode time constants are a function of both 
body pitch and roll damping and rotor tip path plane lag. Control 
over these time constants can thus be exercised through the flap- 
ping lag as well as body aerodynamic damping. 

2. Fuselage 

The fuselage is represented as a virtual flat plate drag 
source having three dimensions. The effective aerodynamic center 
can be located at any position in the body reference frame. It 
would normally be expected to be near the geometric center. 

The fuselage drag model is based on a quadratic 
aerodynamic form originally found in the hydrodynamics text by 
Lamb (Ref erencelO ) and used extensively for airship applications 
by Monk (Reference 11). This form can be easily extended to ac- 
count for fuselage assymetries, lifting effects, and lift 
gradients . 

The simple fuselage aerodynamic form presented here 
provides for drag in forward flight which limits maximum 
airspeed, drag in sideward flight, and rotor downwash impinging 
on the fuselage. All three of these effects are related to power 
losses . 

3. Tail Rotor 

The tail rotor component is modeled in the same manner as 
the main rotor except that no flapping degree of freedom is in- 
cluded. In effect, only Glauert’s equations apply. However 
thrust, induced-velocity, and power effects are correctly 
modeled. Normal directional control is provided through the tail 
rotor collective pitch variation. 

4. Horizontal Tail 

The horizontal tail is assumed to be primarily a lift 
producer, thus only the normal force component is modeled. This 
still provides for computation of drag resulting from induced- 
lift if that is desired. Finally, the effects of aerodynamic 
stall are included. The geometric location of the horizontal 
tail in the rotor flow field is used to obtain the local apparent 
wind component. The location of the horizontal tail provides ef- 
fective static stability and elevator control. 

As with the fuselage aerodynamics, a basic quadratic form 
is used. Two terms model the effects of camber and circulation 
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4 P ne additional term and conditional test 
model the effect of stall. 

5 . Wing 


included 


+ •=4 1 r he component follows the same form as the horizontal 

tain the 3d rela?ed the induced dra * is computed in order to ob 
;fi n .,. the * J related Power- required component which can be 
significant during sustained-g maneuvering. 

o-. Vert ic*I ’P ail 

except T that er f+ Ca i t<>11 13 al f? si,ilar to the horizontal tail 
except that it experiences the flow field produced bp the tail 


C. Definition of Model Equations 

equation"°foJ h fn\ r b! US coraponeI ' t = ~ model are defined, the 

equation*, foi all the components must be expressed in a wav ^ 

ummizes code and the number of parameters The following 

so according to the order of the Lnputer program. f ° ll0Wlng do * s 

1. Main Rotor Thrust and Induced Velocity 

a ciassical C momentum°theory h equation ^uWtri 001 *" ^, ba “ d 

slon scheme which yields" 2“4uick ^genc^" "LTlo'u 

g veT^X® ^ thrUSt ^ indUCad "looiS e.ua^n ^ 



V«ry quickly converging 
Inducad-valoclty loop 
(about 5 iteration*) 


Figure 3 Main Rotor Thrust and Induced-Velocity Block Diagram. 








The recursion relationship is based on breaking the 
thrust- induced velocity loop at the induced- velocity node and 
iterating on a solution for thrust followed by induced-velocity. 
This yields a fast convergence with a fixed number of iterations 
-about 5 is sufficient. 



where 


V W a + <V V “a - b l v a 


W, = W + ^R[0.+7A . ■) 

b r 3 t 'col 4 “twist J 


v 2 = U 2 + V 2 + W (W - 2v. ) 
a a r r i ; 


A = -7T R' 


Once induced velocity for the main rotor has been com- 
puted, one can compute the longitudinal an lateral dihedral 
effects of the main rotor which are, in turn, dependent on in- 
duced velocity: 


db^ /dv = da., /du = 



The main rotor parameters needed for these equations are: 

d , horizontal distance of hub from c. g, 
mr 

h , hub height above the c. g. 

R, rotor radius. 

abcR, product of lift slope, number of blades, chord, and 
radius. 

0 twist ’ e ^^ ec ' t ive blade twist. 

12, main rotor angular rate. 
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^ * Tail Rotor Thrust and Induc©d“V©locity 

Thrust and induced velocity for the tail rotor is computed 
in the same manner as for the main rotor except that no flapping 
effects are included. ^ “* 

The parameters which define the tail rotor effects are: 

d , distance of tail rotor from c. g. 

1 1 * 

h , height of tail rotor above c. g. 


R 


,tr 


3. 


(abcR) , product of lift slope, number of blades, chord 
and radius. 


ft , tail rotor angular rate. 


Fuselage Geometry and Drag 


4 rof lie drag forces are computed for the fuselage in the 
x '. ’ 7 ’ ana 2 axes - These drag forces can constitute a s ; g- 
nif leant portion of the overall power required and thus must be 
computed prior to main rotor torque. The forces are computed at 
the center of pressure located at the point (X.FUS, Y.FUS, Z.FUS) 
relative to the center of gravity. 

Fuselage drag forces are computed using a "quadratic 
aerodynamic form." In this case forces are expressed as a summa- 
tion of terms formed by the product of translational velocity 
components in each axis. The constants in each term are the ef- 
fective flat plate drag. 


w T - w 0 + 

Vj 

local w-velocity 

y ,u » _ P Y f <» 
A »ro 2 A »u 

U.-U, 

drag component 

Y ,u * - P v 

1 wro 2 ' W 

v 0 -v 0 

side-force component 

7 = _P_ 7 ,u * 

C mro 2 ‘-wv 

wr-wr 

downwash component 
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Moments due to the drag forces relative to the center of 
gravity are computed. 

The parameters required for the fuselage are: 

f US 

d , distance of fuselage a. c. from c. g. 
f u s 

h , height of fuselage a. c. from c. g. 

f US 

X , effective flat plate drag in x-axis 
uu 

Y^ S , effective flat plate drag in y-axis 
Z, f , effective flat plate drag in z-axis 


4. Horizontal Tail Geometry and Lift 

The horizontal tail is modeled in terms of a quadratic 
aerodynamic form for airfoils. 


The first 3 tep in computing the lift on the horizontal 
tail is to determine whether the surface is immersed in the rotor 
downwash field. This will influence the local, vertical velocity 
vector . 

The next step is to check for aerodynamic stall by compar- 
ing the force computed above with the maximum achievable at the 
same airspeed. 






L 


stall condition 



k ^ Pitching moment due to the horizontal tail is computed 
center of ° f th * aerodynamic center relative to the 


The parameters required for horizontal tail effects are: 


.ht 


d 1 > distance of horizontal tail from c. g. 


h , height of horizontal tail from c. g. 
Z uu ’ aer °dynamic camber effect 


Z uw . lift slope effect 


Z min ’ sta ^ effect 


5. Wing Geometry and Lift 


, w ?- ng treated in the same manner as the horizontal 

ls checked for exposure to main rotor downwash 

and onen for stall. For the wing, induced drag is computed in 
determine the power loss due to this effect. Lift and 
pitching moment for the wing are also computed. 


wr * w a + V) 

local w-velocity 

c: - i (z:x u 0+ cu 8 vo 

normal force 

> T C.M. U 0 

stall condition 


The power due to the induced drag of the 
based on the product of force and velocity in the 


wing is computed 
x-axis . 


The parameters required for wing effects are: 


jwng di s t ance 0 j> 


wing 


from c. 


g- 
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h , height of wing from c. g. 


wng 

Z , aerodynamic camber effect 

C g . slope effect 

Z™ g , stall effect 
mm 


6. Vertical Tail Geometry and Lift 

The vertical tail is treated the same as the other lifting 
surfaces except that it is assumed out of main rotor downwash. 


. .Vt A . . tr 

V a - V a + Vj local v-velocity 

Y «*ro = T (Y uu u a U a + Y* U a V a ) normal force 
P Vt . 

> “T * r. U 0 stall condition 


The parameters required for vertical tail effects are: 


.vt 


v 

h 


vt 

uu 

( vt 
uv 

,vt 


distance of vertical tail from c. g. 
height of vertical tail from c. g. 
aerodynamic camber effect 

lift slope effect 


3tall effect 




Total Power Required 


Total power ue to 
miscellaneous effects are 
the engine. 


the main rotor, tail rotor, wing, and 
summed giving the total power output by 


Total power required 


pinr + ptr + pfus + pwng + pclimb 


p,mr_ pmr 


,mr 


induced 


jinr 


+ P""- _ . . + p‘ 

profile accessories 


(Note: An estimate of power required for a~ 

scries can be found in Reference 12. 


pmr 

induced ' 


T + v. 

i 


^Profile” " / 2 Cp ° bCK aK[(«R) 2 . 4.6 (U 2 - V 2 ,] 

A d Cl 


C r , bcR 


p"l t _ p, t r 

induced 


T tr 


tr 


,f us 


X, • U 

f us a 


Y- • V 
f us a 


b • ( W - v ) i 
fus a 1 « 


5 wn g_ 


x wng -U 


r,climb : 

1 = m • g • h 


8. Summation of Force and Moment Equations 

i he first order effects of all components are summed in 
three force equations and three moments equations. The force due 
to gravity rotated through theta and phi are also included here: 


X :: - m g sin b ♦ X mr » X fus + X wng 

Y ■ m g sin b cos Y mr + Y tr + Y Vt 

- ~ mg cos b cos 4 Z mr + Z fus + Z ht + z Wng 
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L - 

L mr 

> L fus 

♦ L tr 

M = 

M" 11 ' 

+ M fus + 


N - 

N mr 

+ N tr + 

N vt 

The 

equations of 

motion 


axis accelerations so that they may 
yield body velocities. 


expressed in 
be directly 


terms of body 
integrated to 


9. Integration and Axis Transformation 


As discussed in Reference 13 
numerical integration of states should 
minimise digital effects. 


the algorithm used 
be carefully chosen 


f or 
to 




V 

n+- 


1 


v n + DT(1 - 5 »„ - 0.5 Vl ) 


veloci t 
mat ion . 


Tnese body velocities are then 
ies using a common Euler angle 


converted to earth relative 
direction cosine transf er - 


Finally, the earth velocities are integrated tr -u.. 
position output; 


velocity and 


n+1 


- x t + DT ( 0 . 5 v + 


n 


0 . 5 v . ) 
n- 1 
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10 • Summary of Model Parameters 

A summary of all the parameters included in this model - 

£SitSsr.^K"f. i^s£ r ei 


.n rotor 


F3 . HUB 
WL . HUB 
IS 

E . MR 
I.B 
R. MR 
RPM . MR 
CDO 

A. MR 

B. MR 

C . MR 
TWST . MR 
K1 


Fuselage station of hub 
Watei- line location of hub 

v2!F T& ?? ti 1 lt ° f rotor shaft w.r.t. fuselage 
Effective hinge offset 

Blade flapping inertia 

Radius of main rotor 

RPM of main rotor 

Blade profile drag coefficient 

Blade lift curve slope 

Number of blades 

Blade chord 

Blade twist 

Blade pitch-flap coupling proportion 


Fuselage 


FS . FUG 
WL . FUS 
XUU . FUS 
YVV . FUS 
ZWW. FUS 


Fuselage station of fuselage center of pressure 
Waterline station of fuselage center of pressure 
Aerodynamic quadratic model constant 


Tail rotor 


FS . TR 
WL. TR 
R.TR 
RPM . TR 
A. TR 
SUL . TR 
TWST . TR 


Fuselage station of tail rotor 

Waterline station of tail rotor 

Radius of tail rotor 

RPM of tail rotor 

Blade lift curve slope 

Tail rotor solidity 

Blade twist 


4- Horizontal tail 


FS.HT 
WL . HT 
ZUU.HT 
ZUW. HT 
ZMAX.HT 


fuselage station of horizontal tail 
Waterline station of horizontal tail 

Quadratic max lift coeff of horizontal tail 
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' « 

i 

5. Wing 

Fuselage station of wing 
Waterline station of wing 


Quadratic max lift coeff of wing 
Span 

6. Vertical tail 

Fuselage station of vertical tail 
Waterline station of vertical tail 


Quadratic max lift coeff of vertical tail 


FS . VT 
WL. VT 

YUU. VT 

YUV . VT 
YMAX. VT 


FS . WN 
WL. WN 
ZUIJ.WN 
ZUW . WN 
ZMAX. WN 
B. WN 
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III. Model Matching and Estimation Procedures 

In order to demonstrate model matching and estimation pro- 
cedures, a model of the Bell AH-1S Cobra is developed. The 
actual code for this example version along with a list of symbols 
and a table of associated input parameters are presented in 
Appendices A, B, and C. An example involving the matching of ac- 
tual flight data is presented in Appendix D for the Augusta Model 
109 helicopter. 

The primary sources which are used in the Cobra example 
are the flight manual (Reference 14), a manufacturer’s stability 
and control package (Reference 15), a volume of Jane’s (Reference 
16), and a flight dynamics data report (Reference 17). Other 
useful references include the USAF Stability and Control Datcom 
(Reference 18), the U. S. Army Engineering 

Design Handbook (Reference 19) and the previously cited 
Stepnlewski and Keyes reference. 

In this section the method is described for determining 
the the individual components of the AH- IS and its associated 
parameters. There are 44 total parameters needed for this model. 
22 of these are simple geometrical variables which can be easily 
obtained from scale drawings, from aircraft manuals, or even es- 
timated from a picture of the aircraft. 

A. Mass, Loading, and Geometry Data 

A substantial portion of the data required is either 
directly obtainable geometric data or common mass and loading 
data. 

1 . Geometric Data 

Geometric parameters are easily obtained from aircraft 
drawings or reference literature. Figure 4, taken from the 
flight manual, provides a basis for geometric information. Note 
that positions of all major components are given relative to the 
manufacturer’s reference system (fuselage stations, waterlines, 
and butt lines ) . 

Explicit positions can be obtained for some features such 
as main rotor hub position and tail rotor hub. For airfoils it is 
generally sufficient to estimate and use the positions for one- 
quarter mean aerodynamic chord. The fuselage aerodynamic center 
is less clearly defined and must be estimated depending upon the 
shape, Appendages such as tail boom and landing gear can be con 
sidered i estimating the fuselage aerodynamic center. 




2. Mass and Loading Data 

=*;“E '¥; ;s r 

"sri^rrkE^HL’L^S^; SIsliHSsvi' 

assuming a constant radius ot gyration in each axis resoaled b v 
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TM 66 1620 2M 10 


EXAMPLE 


WANTED 

Om*M**€ APPROXIMATE CENTEX 
OP OAAVmr 

POX KNOWN WEIGHT ANO MOMENT 


OAOSS WEIGHT » 0200 POUNDS 
MOMENT /1 00 - 10100 WCH-POUNOS 

’mcthoo 


move Might prom 1200 pounds 
to approximate nnopowt 
• rrWEIN 10100 ANO 10200 PMJ 
DIAGONAL UNCS 

THIS POINT MOVE DOW* 

TO MEAD 10? 4 ON CENTER OP 
ONAVTTV SCALE 


I 

I 

t 

M 

I 



CENTEX OP GRAVITY PUSIlAG€ STATION 


Representative 
vetght end eg 
chosen from 
flight manual 
operating envelope. 
(90001b, FS 196) 


7»"-IOOM» 1G 


F*eufR 6 t C#nt#* of Gravity L*rr.,u Chxn (Sh*«t 1 of 2 ) 
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Figure 5. Basis for Loading Data. 


r 1 it n n 


ORIGINAL PAGE IS 
OF POOR QUALITY 


Table 4. Basis for Inertial Data 


Condition 


SUMMARY 

TOTAL HELICOPTER MOMENTS OP IKERTIA ABOUT HELICOPTER C.C. 


Weight 

(lbs.) 


ilLi 


Moment of Inertia 
(Slug - Ft 2 ) 

« i •* 7 ^ _• v. 


Principal 

Axi» 

(Nose Down 
Ana la) 


(1) Weight Empty 5571. 4 


(2) Basic 

O) Hog 
(4) Scout 


8673.2 


9501.1 

9296.9 


(5) Most Forward 6606.2 

(6) Most Aft 7476.8 


1990.4 10592.7 .8878.2 

2843.0 13115.6 ' 11264.0 6° 30' 

4002.5 13082.3 11930.4 6* 37' 

3195.3 13233.X U606.71 6° 29' 

2255.5 12462.4 10499.5 7* 19; 

2265.6 11881.0 9903.8 I 4* 18’ 


3 


B. Propulsion Data 

Required propulsion data include power available lor given 
operating conditions. These data can be found in Jane s under 
the appropriate propulsion system manufacturer as illustrated in 
Table 5. The specific information of interest here is the max 
imum continuous power rating for the AVCO Lycoming T53-L-703 gas 
turbine engine. 


Other infornuition needed consists of an approximate break 
down of power, including that due to accessories. Data from the 
Stepniewski and Keyes source are given in Table 6. These data 
will be used to estimate power losses from the computed power re- 
quired by each of the components listed previously. 


The basis for torque (power) available under 
operating conditions is given in Figure 6. (Percent 
assumed equal to percent power for the normal operating 


various 
torque 
rpm 3 


in this case. ) 
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Basis for Propulsion System Data. 


AVCO IVCOMMO GAS TUftWNC ENQWU 


Manufacturer's 

and civil Military 

designation designation 


Type • 


T5313B 



ACTS 

T5317A 

— 

ACFS 

T5311A 

— 

ACFS 

— 

T53-L-13B 

ACFS 

LTC1K-4K 

T53-L-701 

ACFS 

ACFP 

— 

YT55-L-9 

ACFP 

— 

T55-L-7C 

ACFS 

T5508D 

(LTC4B-8D) 


ACFS 

— 

T55-L-1I A t 

ACFS 

LTC4B-12 

— 

ACFS 

ALT 101 

— 

ACFF 

ALF 502R-3 

— 

ACFF 

ALF 502UL-2 

— 

ACFF 


SFC 

T-O Rating M g/J ; % mg/Ns Weight dry Max Length 

tNObst) (Ib/h/hp; less tailpifM dia overSl 

or max kW (hp) 8lb/Mb si) kg (lb) mm (in) mm («) 


Remarks 


1,044 kW (1,400 slip) 
1.119 kW (1.500 slip) 
820 kW (1,100 thp) 
1,044 kW n,400 thpl 


1,157 kW (1,550 shp) 
1.082 ekW (1,451 ehp) 
1,887 ekW (2.529 ehp) 
2,125 kW (2,850 shp) 
2,186 kW (2,950 shp) 
flat-rated to 
1.678 kW (2.250 shp) 
2,796 kW (3,750 shp) 
3.430 kW (4.600 thp) 

7 2 kN (1,620 lb) 

29 8 kN (6,700 lb) 

33 4 kN (7,500 lb) 


98 (0*58) 

99-7 (0-59) 

1 15 * 0 - 68 ) 

98 (0-58) 

98 7 (0-584 r 


1014 (0 60) 
102 7 (0 608) 
101 4 (0 60) 
100-1 (0*592) 


245 (540) 
256 (564) 
225 (496) 

W 

312 (688) 
363 (799) 
267 (590) 
274 (605) 


584 (23) 1,209 (47 6) Powers Bell 205 A 

584 (23) 1,209 (47 6) Based on T5319A 

584 (23) 1,209 (47 6) BeU 204B 

_384 j23) U09 (47 6) Advanced UH-1H, A H-lG 


1*4 


106 0 (0*628) 

89 6(0 53) 322 (710) 

86-2 (0-51) 329(725) 

110-19(10 36) 156 (343) 
811-64 ( 10 41 1)563 (U45) 
812 1 ( 80 428) 590(1,298) 


1 * 4 , 

584 (23) 
615 (24 2) 
613 (24 2) 
610 (24) 


615 (24 2) 
615 (24 2) 
584 (23) 
1.059 (41 7) 
1.059 (41 7) 




1,483 (58 4 
1,580(62-2) 
1.) 18 (44) 
1,118(44) 


1,181 (46-5) 
1,118(44) 
890 (35) 
1.443 (56 8) 
1,487 (58 56) 


fcell AH- iOTaR^ S TtSW/Cobrt ! 
"BeU XV- 1 J — 1 


Grumman OV-1 D 
Piper Enforcer 
Boeing CH-47B, Bell 214A 
Bell 214A. 214B 


Boeing CH-47 
Improved T55-L-11A 
NASA OCGAT 
BAe 146 

Canadaw CL -600 Challenge* 


t A^L,’ i^T*55 P ! U ‘i Ta 1 ’ n ’if"* P ,u ‘ eenlribii»l, bee-turbioc propeller; ACFF - u.j plus cen,nfu»»l. bee-turbue f, 

TAppket to TJS-L-IIA. C , D. E end 712 , thotc deufnMed * * bavuf 2H ain coalinfeoey rating of 3,357 kW (4,300 thp). 


Table 6. Assumed Breakdown of Power Absorbtion. 


% 

Total Power 
in Hover 

% Total Power 
Max Forward 

Main rotor induced power 

65 

15 

Main rotor profile power 

15 

50 

Fuselage parasite power 

5 

05 

Tail rotor total power 

10 

5 

Misc . and accessories 

5 

5 


( NOTE : 


Power losses due to wing stall should also be con 
sidered where the effect is suspected to be 
significant. It will be neglected in this ex- 
ample . ) 
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Rotor Data 


Rotor system characteristics consist of geometric, 
aerodynamic , and operating condition features. Most of the 
geometric data including sise and number of blades and hub center 
a*e easily tound in flight manuals. Operating conditions, namely 
the normal operating rpro, are likewise obtained. 


The main aerodynamic 
tion lift curve slope and 
accepted values of 5.7 
starting points 


parameters include the effective see 
profile drag coefficient. Commonl 
and .006, respectively, are sufficien 


The most crucial rotor parameters. however, are those 
relating to the effective flapping stiffness or hinge offset. 
These data are generally found only in manufacturers design 
reports. Of course in the case of a simple teetering rotor the 
effective hinge offset is zero. Articulated rotor designs are 
also fairly easy to represent as long as the geometric hinge of- 
fset is known. The most difficult variety to model is the 
hingeless rotor since both an effective hinge offset and flapping 
spring must be determined. 


Useful auxiliary information for modeling the rotor system 
is response data which provides direct indication of the unaug- 
mented pitch and roll damping. 


D. Aerodynamic, Features 

Aside from the rotor system aerodynamics, parameters must, 
be estimated for the airfoil and fuselage components. The tech- 
niques for doing so are common and require little effort. If 
manufacturer’s stability and control data are available these 
calculations are trivial. Otherwise, one can refer to estimation 
nandbooKs such as the USAF DATCOM (Reference 18) . 


b- r 
a re 
ing 


Airfoil Lift parameters involve three main features: cam 

and^ incidence circulation, lift, and stall. The first two 
highly dependent upon geometry and the third on maximum lift 
perf ormance . 


Relationships which are needed for setting parameters in 
vi’ive the quadratic aerodynamic parameters and the more common 
non dimensional aerodynamic coefficients. These are given below 
for use in the estimation procedures described in Figure 7. 
The equations are for the horizontal tail, but the other airfoil 
surface:; are similar. 


Kst.imates typical for airfoils: 
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2 “ v = - S ht d? 1 '’ ; C T is set by both camber and incidence 


,ht 

J uu 


'L ’ L 
o o 


of the airfoil . 


,ht 


Z“ w = -S ht C ^ ; note that C T x 
uw L a *->«< 


2 jJK 


M+2 


Z““ = - 5 ^ ; typical values are 1.5 to 3 depending 


,ht 

J min 


max 


upon aspect ratio. 


Similarly, fuselage drag estimates can be made for each of 
the three axes using available drag data. 


Estimates typical for fuselage drag: 


v fus_ r fus 

A — O 

C D 



UU 



, n f us . 

where b is 

the 

projected frontal area 


and C D can 

be 

estimated using numerous 

text boo K 


tabulations of 3 -dimensional drag. This will 
vary for each axis. 




<j 

• Ml 



WING: 

span= 10.75' (based on i = 14\ 
chord = 3.0' *1.2, 

area = 32.25 ft 2 assume C L » 2) 
aspect ratio = 3 m “ 

C, * hi & = s 
At * 2 


HORIZONTAL TAI L 
span = T 
chord = 2 25‘ 
area = 16 ft 2 
aspect ratio = 3 


TMM.1M0-2M.10 


VERTICAL TAIL: 
span = 5' 
chord = 3 3’ 
area = 17 ft 2 
aspect ratio = 1 




(assume C = 


— I«f?| — 


■T'i 




. si rr i m _«**»«**« U»cth 

J *OT0*S TUAN f»c 




mm at*!’ *** FO * ** ,0 « f i*ct«o 


FUSELAGE 
assumed drag 
coefficients 
front, 0 2 
side, 2 0 
plan, o 7 


"**•*-* MndMOtaMiMi 


(Two D 


Figure 7 


Basia for Initial Estimates of Aerodynamic Parameters. 


E. Hover Performance 


The parameters listed above provide a starting point for 
the math model. Additional flight manual and available flight 
data will serve to make refinements in model response and perfor- 
mance characteristics. 

The first adjustment of model parameters can be made based 
on the flight manual hover performance as shown in Figure 8. 
Here the percent maximum torque is given for a specific hover 
condition . 

The factors which can be adjusted to achieve a good match 
are the power losses due to accessories, downwash on the fuselage 
and horizontal airfoils, or main rotor induced velocity factor 
(if included). 
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TM U1I20 2M10 


Torque req'd 
for OGE hover. 
< 1 232 hp) 


HOVER 

1 CONFIGURATIONS 100% PPM 
UK VC l SURFACE CALM WIND 


MOVER 
AM IS 
TltlTOJ 


EXAMPLE 

WANTED 

TOAOUi RE OUlAE D TO HQVE A 

ft MOWN 

pressure AtTnuot * iioooKrr 

FAT • 0 C 

GROSS WEGmT = SSOOLA 
DCSIAEO SKID HEIGH? * 2 FEET 

MET MOO 


ENTER PRESSURE A*. TIT UOC KEA* 

MOVE RIGHT TO I AT 

MOVE DOWN TO GROSS WEIGHT 

MOVE LIFT 10 SK*0 HUGH* 

MOVE DOWN READ TOROuE R£0WA*9 
TO MOVER * 72 SVJ 



TORQUE -%Q 


DATA AASIS Of RIVED FROM F UGH 7 T|$T 


I 10 IS 20 

OCNSITV AlTlTUOt 1000 FEET 


8 

■ 


Flours 7-§ Hov«r chart (Sh««t 2 of 2) 


7 17 


Figure 9. Basis for Hover Power Required. 
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Forward Flight Data 


f to thls point model adjustments have centered on the 
moiii rotor system since body drag has been low due to the hove,' 
oor.di tion. With the consideration of forward flight the fuselage 

iow plays a major role in limiting maximum speed and climb re*-- 
* o i m .4 Hue . ■ ^ ^ 


. ... . The . m ^ in set n of data useful for adjusting fuselage drag 
a,e given in figure 9 from the flight manual. Note that the 

fnThi 7 information is the torque required as a function of 
flight condition and loading. The two main features on this riot 
t ’_ ' ' ie maxim um speed at continuous operating torque and the 
toique and speed for level flight at minimum power 


Additional information is given in Figure 
maximum rate of climb corresponding to an increase in 


10 with 
torque . 


1 hrr 


Finally in Figure 11 data are driven for the 
and minimum rate of descent. These are useful for 
effective full-down collective pitch stop. 


max imum <? ! \ de 
sett i ru r r .he 
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TM M-1620-236-10 


Torqu* for imx 

level speed. 

( 133 ki # ease) 


Torque and speed 
for level flight at 
min power. 
<64kt *46%) 


CRUISE 

PRESSURE ALTITUDE — SEA LEVEL TO 2000 FEET 
100% PPM, CLEAN CONFIGURATION JP-4 FUEL 


2000 FEET 

PRESSURE ALTITUDE 


CAUISC 
AH IS 
TS3 l 703 



liMLHr Ft9ur» 7 7 CruiM ChAd (SNmi 14 of 23> 


Figure 10 . Basis for Forward Flight Speeds and Power Required. 
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Aft a final note, the process of tuning mode] parameters 
..hould not be done without careful consideration of all secondary 
effects. The best policy is to avoid making anything other than 
simple direct first-principles corrections. There is substantia] 
redundancy in some of the data shown here, and it is not possible 
to achieve perfect matches in all respects. One needs to exer- 
cise judgment in the degree of accuracy required as a function of 
the model application. 


I v Checkout Procedures 

A. General 


thoroughness of simulator 1 ^©^? r 1 Complexit y can hamper the 
checking However the moWei put r Program implementation and 
with reasonable Effort* 6 ?m* f r “* nt ? d h ? re oan *>« fully checked 

constants and degrees 'of freed au ® !° V** 3l " a11 number of model 
degrees of freedom, and minima] program branching. 


Th.a 


recommended checkin/? procedure involves 


merits 


the foliowin/z' 


a. 


• Use of an independent operating program. 

• Verification of trim points. 

• Verification of state transitions through n step 

• Overlay of time histories. 

• Identification of dominant response modes. 


oome of the sp ^i-^ncr 
t, - ij a . , 1 si:eps an-s reuunaent 

to build confidence in the correctness 

piemen tation at only minimal added 

brief discussion of each element. 


^ w u 


* a. v j 




the math model im' 
cost. The following is a 


li. 


Di 


scuss ion 


of Checkout Procedure Element: 


1 . 


ndt 


^pendent Operating Program 


comp 1 i shed us i ng"^ independent *,!! , model checkout should U. as- 
Furthermore. not only Tivoli and check sou, -e. 
also an independent computer *" independent Program be used b.-t 


This 

model version 


m on h a small JT? enableS the us * r to develop a math 
sets of check cases win ? ^ microcomputer and run complete 
outer facilities. " n advance of using the simulator com 

math model consisted S ^Compa^JflG i®®!.* 0 develop and run thi,. 
working memory running 'Micro-eft 286 desktop co "‘Puter with 640K 
was used aUhZughTK^ interpreter mode 

Permits a highly efficient 'interaction interpret, 
«.e\ eloper and the computer system. 




between the mode 


2. Trim Point Verification 

A check of static trim points gives an initial indication 
of correct model implementation. The full operating envelope can 
be covered with just a few cases and possible di screpenci.es iso- 
lated to airspeed. vertical velocity, or controls. A cursorv 
check of suspected parameters or component equations can usually 
lead to simple corrections. Trim solutions should be correct 
prior to proceding to the next item. 

A sample of the trim solution printout is given in Figure 
12. This same format is displayed during the trimming process 
so that one can observe whether there are difficulties in iterat- 
ing on a solution. 


"Rift CALCULATIONS HEFHEL2: FULL UTILITY VERSION 

CONFIGURATION: 102 AH-1S 
05-25-1987 16:09:55 


Foot = 1.05E+00 

DC 

= 

15.7 

fidot = -1.44E-01 

al 

= 

1.3 

Rdot = 3.28E-03 

bl 

r 

-2.1 

Udot = -4.B8E-02 

DTR 

= 

1.02E+01 

Vact = -3.87E-02 

Theta 

= 

-1.3 

Ndot = 1.41E-03 

Phi 


-1.02E+00 

aldot= -6.01E-02 

Bl 

- 

-1.30E+00 

fcldot= 9.80E-02 

Al 

= 

-2.05E+00 

Q = 1 . 34E+04 

HP 

= 

973 

Vi = 35. B 

Thrust 

; = 

9256 

Vi.tr= 47.9 

T.tr 

= 

618 

VBI1!= O.OOE+OO 

Xdot 


O.OOE+OO 

VB (2>= 0.00E+00 

Hdot 

= 

0.0 

VBI3)= O.OOE+OO 

Gama 

= 

V.0OE+00 


VT a 0.0 

Hit 10) to Freeze tna anvtue 

Triaied: Hit PRTSC to aale hard copy 
Hit RETURN to continue 


Figure 13. Sample of Trim Point. Printout. 




3. State Transition Verification 

Given that static solutions are valid, the dvr, 
espouse characteristics should be examined next. Correct op 
•ion is indicated by tracking several discrete state van 
transitions and comparing with independently obtained c 
frrii-F J h , s ma< ? e fea sible by restricting the number of 
"•x ample onlT^ t®"- ] evels of numerical integration, 
needed ^to 2t y ^? b 1 ^ "transitions for each control variable 

e to exoitt; each term in the model equations. 

In order to thoroughly check state transitions a t 

"ta? '\n l ' e V° mmenr!ed - This is accomplished bv duplicating 

that -f +h ^ XO ? . printout ’ forDiat of the checkout computer 
that Of the simulator computer. The original checks ,n 

printed on transparencies then directly overlaid with 
simulator printout. 1 " 


h o c k 
d e 
For 


a t, ] e 
the 
wi t.h 
he 
the 


Examples of the stat.e transition checks 


are given in Table 
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Time History Overlays 


Y 


y W 


4 . 


In theory the combination of static and state transition 
checks should be sufficient to demonstrate agreement with the in- 
dependent model implementation. Howeverg additional confidence 
is gained by selecting several time history cases to overlay. 
These can be supplemented by checking dominant response modes 
based on transfer function solutions from the original independ- 
ent check model . 

Useful time histories to consider are angular rates for 
both on- and off -axes for a given control input. This checks 
both the dominant response modes and the amount of off-axis cross 
coupling. Examples are shown in Figure 13 corresponding to tn- 
previous check information. 



t 

i 




Of- i'Wn 


• ' ' J . r S 

•«u run 




I 
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I ,oi* l aiva HDiid 


I ,01 ■ I n v a m 



I ,o !♦ I 3iva noa ' " I ,oi • I iHd 


Figure 13. 


Examples of Time Histories to be Used for Overlays 


TIME 



JF tn -r ‘ •' & 

F ' 0c *< quality 


Dominant Response Identification 


It is also useful to supplement the above checks with a 
comparison of identified dominant response features from t he 
simulator computer with those features observed or computed i r .,m 
the independent checkout version. This is particularly important 
for handling qualities investigations. 


Dominant modes are examined bv exciting an axis with tt.» 
corresponding direct control and scaling the appropriate first- 
cr second-order response features from the respective motion 
traces . The on- axis traces presented earlier in Figure 13 serve 
t iis purpose for extracting short-term pitch response informa - 


t i o n . 


- 53 - 



V Model Extensions and Refinements 

The example which has been presented above can be modified 
in a number of ways in order to address specific simulation 
needs. The above math model can be either simplified or made 
more sophisticated. The following is a discussion of some pos- 
sible extensions and refinements. 

A. Flight Control System 

There is no flight control system included in the above 
model other than conventional aerodynamic interfaces such as 
cyclic, collective, and tail rotor controls. Addition of a 
flight control system requires definition of relationships be- 
tween the cockpit manipulator and the above aerodynamic controls 
plus any stability and control augmentation systems. 

As with the basic airframe math model, definition of 
flight controls can be done with a wide range of eomputati uiiai 
eomulexi t.y . However the same considerations can he applied in 
•order to match the level of complexity with user utility. The 
main question is to what degree can the simulator pilot observe 
or be influenced bv math model intricacies. 

B. Engine Governor 

This aspect of the helicopter math model can be important 
for tasks involving maneuvering or aggressive control of collec- 
tive pitch. 

The above math model, is designed to accomodate an engine- 
governor system since rotor speed is explicit in the equations. 
It is necessary only to add appropriate engine governor equations 
of motion prior to computation of the main rotor thrust. 

In general, only a second-order engine governor response 
is required in order to handle the effective spring-mass-damper 
action of the main rotor combined with the propulsion system and 
governor control laws. An adequate model is described in 
Reference 20. 

C. Ground Effect 

The modeling of ground effect can he important for tasks 
involving hover under marginal performance conditions. Again, 
the computational complexity of such models can varv widelv. 

f t, is recommended that, as a first, out. ground effect be 
iVii.de 1 , ed as an induced-velocity efficiency factor which primarily 
if tecta the thrust, arid power required to hover. This efficiency 
factor can be adequately modeled as an exponential function of 
altitude. The exponential scale height and magnitude is easily 
quantified from the flight, manual hover performance shown earlier 
in Figure 8. 
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nc'i& L fS 


. dynamic Inflow 


° f p 00« QUALITY 


oort.ant/to" model^the ^eo??ie r ur?n e t .; PPl i C 2 tlons 11 ■** in.- 
Pit.h chan*.. This u ?ypIcaUy a 


the < 

l)r- t'j ( 


moments 
R*?f t*renot 


T}iis effect can be modeled bv setting a fi r = t-, yH & , , 
y* ou I at ion of thrust and induced ve loci tv i^f - : J * £ 

resulted for guidance in setting values ' ‘ «**«*K:e ‘ A ' V)n 

also 1 - 


rj 


:an 

n n 


, ** ~T V"'”" vajuca - Other forces 

be affected by dynamic inflow as described 


■i r i * 1 
i n 


E. 


Higher-Order Flapping, Coning, and Lead-Lag Dynamics 


exaainin** h f ii*ht <r C ontro 1 * V avat d "“? lC * mav b * ° f 

effects . ^ However there's ' 

tional ability of the* + ^ out^iae the eamonta- 

system. Thus is c^ucill fo? ^LlodefeJO 1 sto , rted ^ ‘he motion 
requirements relative to capabilities 1 ^ c ‘ naly2e computational 
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APPENDIX a 

BASIC PROGRAM LISTING OF MATH MODEL 


‘ 480 
: 490 
*500 
1510 
1520 
1530 

J 540 
* 

» . j . 

? 5 7 0 

15S0 

♦ rc'i 


lb2‘) 


. sal’ 
1630 

, - V(: 
! '■ 0 
1710 
; 720 
173- 
:740 
1750 
17i0 
17 7 C 
1 7EC 

i “n; 

M'jf, 

' 510 
:S2': 

• J30 

1250 

*260 

1£’(> 




Preliiinarv Calculations 

V7.IN.rps = ABS ( X DOT/COS ( 6AHNA . RAD ) j 
VT - VT.IM.FPS/FPS. PER. KNOT 
V7.1N.P?S.?3UARED * VT.JN.FPS - ? 

*1 = WT/5RAV 

0K5A.MR = RPK. MR*2«PI/60 
0SEEA.TR = RPM.TR»2*PI/60 
V. T !P=R.NR*0MEGA.MR 
C F:.NR * CD0*R. MR#B, MR#C. RR 
■ r F:.^ = CBC»R.TR«B.TR«C.TR 
TP. LOSS = 90 

fr “ *■«»* 

PRESS. RATIO = TENP.RATIO-TEMP.EXP 

JENS. RATIO = PRESS. RATIu/TEUP. RATIO 

qo r OEMS. RfiTI0*RH0. SEA. LEVEL : R2-R0/2 

6AM.0M.16 = R0*A. MR»C. HRfR, HR A 4/I.B *0HE6A.NR/16#U + B/3*E.MR/R HR) 

- 5 ! • IOMEGA. MR*E.MR/R. MR / 6AM OH 16 ) ♦ n.- t\ ‘" R 

i.sis*i=s.».t.KR.c.m/R.»/pi .. lcl! " ! 

Tpf «HKt 

: ' TPP oitchup with speed 


-1. HUB 
A. "OS 
t.NN 
H.HT 
'’.VT 
H.tp 


. WL . HUB-WL . CS) /] 2 
XWL.FUS-W1_.CG) /!2 
{ WL.WM -WL.CG)/12 
(ML. IP -ML. 05); 12 
'ML.VT -ML.C61/12 
(NL.TR -ML.CO/12 


D.HUB * (FS.HUB-FS.CG)/12 
0.FUS = (FS.FUS-FS.CG)/12 
D.NN : 'FS.NN -FS.C6I/12 
IPS.HT -FS.CBI/12 
(FS.VT -FS.CB)/12 
(6S.TR -FS.C8I/I2 


0.HT - 
D. VT : 
D.TR = 


hub re eg 
Fuselage re eg 
' wing re eg 
horizontal tail re eg 
vertical fin re eg 
tail rotor re eg 


A-l 


ocT'JPN 


I*** N ■ 




EViwHirE: Dynamics subroutine 

1 ... ( .»hhhh»h«#h > hhh«mhhh»h*»ihihhm»hhhhhhhh 

4vtm«4«*»4m Prelieinary calculations *»******«##**» 

C4 - CCSiXEWI : S4 = SlflliXE (4) ! evaluate Euler angle trio fns 
»- = CCS (XE (5)1 i S5 * SIN(XE(5)! 


VA; 11=73(11 
V'A (21=73(2) 
VA!3!=VE(3) 
VA;4?=VB(4) 
VA!5>=78(5) 


-CVStli «C5) evaluate relative aireass velocities 
-(7S(2)*C6-7E(li#36) 

-(VGI3)«C5+Vfi(l!*S51 

-VB(4! 


•'T,i - St-f: ! 7 A ! 1 ! » V A 1 ! + V A ( 2 ) ♦ V A ( 2 ) ♦ V A ! 3 ) * V A ( 3 ) ) 

hohhhhhjh Rotor tip cath plane dynamics >ihh«mhh 


A.;lK= GV(8)-DC(2)+KC*SV(7)+DB1DV#VA(2) 
5-7(1= SV(7)+DC(3)-KC»GV(B)+DAlDU*VA(lf 


bl - A1 + e.al + dbl/dv .V 
al + Bl - e.bl + dal/du .U 


GA(7)= - 1TB#B.SUH - 1TB2.0H4A.SUH - VAi5) a!. dot 

SR(8S= - ITB4A.3UH + 1TB2.QH4B.SUM - VA‘4> bt.dot 


Sv*7)=6V(7? + ST*(A2*6I1I7) + 82»AP(7)) 
5V(e»=5V(e; + ST«IA2«6R(B) 4 62#AP(8») 

A p 7)=GR(7> : AF!8)*SR(8! 


al updated 
bl updated 

save past values 


♦4*» >444444 Main Rotor thrust and induced velocity ** 4 t 44*#***»4 

*1 " ^ - IS),VA{1! - GV(BI4VA<2) ;• e-axis vel re rotor plane 

*1 * m *2; tiA. HR4R.HR* (DC ( 1 1 * .75#THST. HR): 'e-axis ve! re blade 

‘ “ ,l ' ■' w : iterative solution of thrust and induced vel 

TH\"JST . HR= («-V! . HR) *DHESA. HP#R. HR»RHO#A,HR*B. HR*C. HR*R. HR/4 
7HAT.2=VA'.l ) A 2 4 VA(2) A 2 ♦ WRf (NR-2*VI.NR) 

v !. HR. i =SQR ( (VHAT, 2/2)4 (VHAT. 2/2) * (THRUST. HR/2/ (RH0*PI*R.HR A 2)) A 2) - VHAT 2 >"> 

v:.HR=SBPiABS<7!.HR.2» .a, n rotor induced velocity 

^ t / T J 

„M**M444VM44*t4*l#m»#, FuselaO* *4*t*»4*l*l, *4, 

= VA!3) - VI. MR .• mdude rotor dounwash or. fuselage 
D.r*=. Va U - / ( -W. FUS) * IH. HUB-H.FUS) )-( D.FUS-D.HDB ):' cos of dounuash on fus 

> - FU£ = R; * : ™.FUS » A8S(VA(1) l * VA(1) draq force 

’ P * m,FUS * ABS(VA!2)) * VAf2) side-force 

-■'-u'S = R2 * 2HH.FUS * ABS(WA.FUS) » HA. FUS heave force 

l.F'JS * T.FUS*H.FUS 


ORIGINAL PAG5 is 
OP POOR QUALITY 


ORIGINAL PAGE IS 
OF POOR QUALITY 


> 



" : 40 


-:i0 




4 r 0 




4150 

4: t : 

;:?0 

-290 
4300 
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H.FD2 = Z . FU5*D. FW - X.FUS«H.FUS 

MKDLCED.HR = THRUST.HR » VI. HR 
P.CLIN? = WT*HDDT 

'-.PARASrZ = - LFUS»VA(1) - Y.FUS«VAI2) - Z.FUS«MA.FUS 
F. PROFILE. HR = F2MFR.HR/4)«0HE3A.HR*R.HR§(0HEGA.HR A 2*R.HR A 2 ♦ 

4. :* iVA( ! ) *VA ( 1 ) +VA(2) *VA ( 2) ) ) 

POO.HR = P. INDUCED.HR + P.CLIHB + P. PARASITE ♦ P. PROFILE. HR 
P5XER.FDT0R.HR * P. INDUCED. HR + P.PROFILE.HR 
MEMOS - P.FARASITE 
'DR5JE.HR = F9HER.hr/0HEBA.HR 

CofliCLite main rotor Force and aoaent coapcnents. 

4. HR -THRUST.HR * (6VI7HS) 

V. M R - THRUST.HR ♦ SV(B! 

:.hs = -thf.ust.hr 

L.HR - Y.HR*H.HUB + DL.DBt#6V(B> + DL.DAHKBV(7HDC(3) - K1«6V(8)) 

-•HP - 2 .!1R*D. HUB - I.HR*H.HUB + 0L.DBUBVI7) ♦0l.DAl#l-GV(B>+DC(2) * K1#6V<7! > 
N.HR = T0R5UE.HR 

**«mm*tt«* Tail Rotor thrust and induced velocity ***♦«♦*♦*#**« 

rR.TR = - { VA f 2) - VA<6)»D.TR + VA(4!*H.TR; velocity relative to rotor plane 
VP. TP = VR.TR +2 /3«0HE6A. TR*R . TR* (DC < 4) +THST . TR*. 75! velocity relative to blade 

FOR 1=1 TO 5 iterate on thrust and induced velocity 

THRUST. TRs(VE. TR-VI . TR) fOHESA. TRtR. TR*RH0*A. TR#S0L.TR*PI*R.TR*R. TR/ 4 
VHAT . 2* ( VA (3) +VA ( 5) *D. TR! A 2 + VA(l! A 2 + VR .TR*IVP.TR-2»VI.TR) 

VI . TR. 2=SBR i (VHAT. 2/2) ♦ (VHAT. 2/2) + (THRUST . TR/2/ (RHQ*PI*R.TR A 2) ) A 2i - VHAT. 2/2 
VI,TR=S5R(ABS(VI.TR.2! ! 

HEKT : 


PQXER.TR = THRUST. TR*VI.TR 
•MR = THRUST.TR 
..TR = Y.TR*H.TR 
N. 'R * -v.TRaD.TR 

»»*<*********»*«+#** Horizontal tail *****f*«*m»*»**mm***m* 

•.DX=! VAIi:MVI.HF:-VAi3)*(H.HUB-H.HT) )-( D.HT-D.HUB-R.HR ):' dnash iapinoes on tail? 
EPE.HTs.5»(:*SSN(D.CH)) uni (ora doanaash Field 

IF D.D»')0 AND D. DIKR.HR; THEN EPS.HT* 2*(H).0H/fi.HR) ELSE EPS. HT=0 :* trianqlr dnash 

= VAC! - EPS.H T *VI.HR + D.HT»VA(5) : local z-vel at h.t. 

'.'TT.HT=SQR(VA(1!*VA( l'+VAi2)*VAI2)*NA.HT*NA.HT) 

- • HT-R2*!ZUU.HT*ABS(VA(ij )«VA!1! + ZUN.HT*ABS(VA(1))*NA.HT! circulation lilt on h.t. 
IF AB5;«A.HTi .3«ABS(VA(l)i THEN 2 . HT=R2* ZHAX . HT*AB5 ( VTA. KT > *WA. HT surface stalled’ 

H.HT = Z , HT*D. oitchino aoaent 
omH******************* Hinq 



1 


i 

1 

r 

* 

\ 

i 


i 




WA.WN = VA * 3 ) -VI .MR local i-vel at »ing 
VTA. WN "SfiR ( VA ( 1 ) *VA ( 1 1 +NA. NN*WA. UN) 

2 . W'N=R2* ( 2u‘U. HN*VA ! 1 ) *VA ( ! ! ♦ ZUN.NN*VA( 1HNA.NN) nortal force 

*.M a R2/F I y VTA. WW/ VTA. WN* < 2UU. «N*VA ( 1 1 # VA (1 i +2UM. HN#VA 1 1 ) #WA. MNi A 2 : ' inouced drag 

- r A35!WA.WNi>.3»ABSlVAIl!) THEN Z.NN=R2*ZNAX.NN*ABSIVTA.WN)*NA.WN surface stalled? 

POK'i.Wi = ABSiX.WN*VA(!)i 

POWER = POWER. MP + PONER.TR ♦ POWER. HN ♦ HM0SS*550 

H*u#**mmm**H*« Vertical tail 

Ccmpute aerodynamic forces on vertical tail 

VA. VT-73 12; -*VI , TR-D. VT*VA ik) 

VTA. VT=5GR ! VA ( 1 ) *VA (1 ) + VA . VT*VA . VT 1 
. v'T* Q 2* T YU'J. vT»ABS ! VA ( 1 ) 1 * VA 1 1 ) ♦ YUV.VI#ABS(VA(1))#VA.VT) 

;F ABE VA . VI ! .• . 3*A8S ( VA ! ! ) 1 THEN Y. VT*R2«YHAX. VT*A8S ( VTA. VT) *VA. VT surface stalled? 

L.' ;T = y.vt«h.vt 

N, VT = -Y.VT*D.VT 

«*v*f»»m**m* General force equations hkhiihihihihhihh 


. ERfiV * -H*BRAVtS5 
. SF.'AV = K*GRAViS4«C5 
. 3RAV * K«GRAV*C5*C4 


: gravity forces 


• gravity 1 H.R. 1 FUS. 


T.R ! H.T. ! WINS ! V.T. 


conoonent 


.1; = X . BRA V + X.MR * X.FL'S * X.NN 

r( 2) = Y.BRAV + Y . HR + Y.RJ3 + Y.TR * y. VT 

F!’: = 2.GRAV + Z.MR * 2. FUG + + 2.HT ♦ Z.NN 

F!4; = * U. MR ♦ L.FUS + L.TR + L. VT 

F (5: - ■* H.HR + M. FUS + H. HT 

Flat + N. MR r N.TR + N.VT 

1 1 1 : t i I 

-(?; * GR<7) /ITP 

C 'S) = SRiS’/ITB 


X-force 

Y-force 

2-force 

L-»e«ent 

N-ioiert 

N-»o«ent 


pitch flao 
roll flao 


IF CHECK*® THEN GOSUB 7790 : ' fill force co.ponent array 


Body Accelerations 


ABC! =• - i VB (5 1 *VE i 3) - VB ( 6 ! * VB i 2) ) ♦ F'll/f! 

ABC = iV6'4WB!3)-VB<mVB<6)) ♦ FI21/N 

ABI3 1 = !VBU)*VBi5!-V8<41*VB(2>! + FC!/H 

A3t*> = F(4!.TX 

ABC; = F'E. 'IY - VB!4l«VB(6!»(IX-I2)/lY ♦ (VB(4)»VBI4>-VB(4)*VB(4"«I)(2/IY 
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52 60 
z::o 
5280 
52*40 
5300 
57! C 


j -.‘ i 

:?c 


c • "V. 


SBli) * F 61/12 * IK2*fl8(4)/lZ 
:nlrirat£ Body A:ce! erat l ons 
F0? !?. = ! T 0 £ 

VB !I; - VP (IX) + ST f (fll ♦ AB(IX! + B1 * AF'(IX)) 

« AB(JX) : REH SAVE ACCEL PAST VALUES 

ix 

Trar.=fora to earth (A/C re! to deck) velocities 

VE.l! - (VB(1 ) # C5 + VB (T) * S5) ♦ C4 * COS ( XE (6> ) 
vE (2) = VB (2! *CQS ( XE (6 1 ■ +VB ( 1 ! * SIN (XE(6) > 

VE!3! * (VB(!; # S5 - VB(3! # C5 ) # C4 

V‘E'4; = VB(4) ♦ 108(5) » S4 + VB(6) * C4) « TAN(XE(51) 

VE.5) = VETO * C4 - VB(6! * S4 
VE'S) = 1 VB!6i * C4 1 VB(5» * S4) / C5 

Integrate earth (A/C relative to deck! velocities 

rC 1 ' 1 IX * ! TO i 

SE'IJ* - XE(IX) ♦ ST * <A2 * VE(IX) + 82 * VP(IX>) 

VFCr.) * VE'IX) : REH SAVE VEL PAST VALUES 
h£< T IT. 

T I1E-T!HE+ST 

IF CHECK =1 THEN IF CHECK. LQ0P< CHECK. LOOP. MAX THEN GOTO 3520 
! c CHECK* 1 THEN GOSUB 8890 

RETURN 
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APPENDIX B 

DEFINITION OF PROGRAM SYMBOLS 


r 






AB(i) 


AB( 6 ) 
AB( 7 ) 

AB(8) 

AP( i ) 

A. SIGMA 


DAI DU 

DB1DV 

DC(i) 
DC( 1 ) 
DC( 2 ) 


Numerical integration constant (Adams- two = 1.5). 
Numerical integration constant (trapezoidal = 0.5) 
Body-axis acceleration vector. 


AB( 1 ) Body x-axis acceleration (U) component (ft/sec 2 ). 

AB(2) Body y-axis acceleration (V) component (ft/sec 2 ). 

AB( 3 ) Body z-axis acceleration (W) component (ft/sec 2 ). 

AB(4) Body roll axis acceleration (P) component (ft/sec 2 ). 

AB(5) Body pitch axis acceleration (Q) component (ft/sec 2 ). 


Body yaw axis acceleration (R) component (ft/sec 2 ). 
Lateral tip-path-plane angular rate (b, ) (rad). 

Longitudinal tip-path-plane angular rate (aj) (rad). 

Past value of AB(i). 

Product of lift-curve-slope and solidity. 

Numerical integration constant, 1-A1. 

Numerical integration constant, 1-A2. 

Cos[XE(4)] or Cos of roll Euler angle. 

Cos[XE(5) ] or Cos of pitch Euler angle. 

Cos[XE(6)] or Cos of yaw Euler angle. 

Thrust coefficient. 

Partial of longitudinal flapping to forward velocity 
( rad/ft/sec) . 

Partial of lateral flapping to side velocity 
( rad/ft/sec) . 

Control vector. 

Main rotor collective pitch angle (rad). 

Pitch control, B^ , (rad). 
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DC(3) Roll control, Aj, (rad). 

DC(4) Tail rotor collective pitch angle (rad). 
DENS. RATIO Density ratio. 

D.FUS Fuselaae horizontal T*r \ a 4 + 4 


D.FW 

D.HT 


D. HUB 


F(i) 

F(l) 

F(2) 

F(3) 

F(4) 

F(5) 

F(6) 

F(7) 

F(8) 


Fuselage horizontal position of aerodynamic center 
(ft) . 

Position of downwash on fuselage (ft). 

Horizontal tail aerodynamic center relative to c. g 
(ft). 

Hub horizontal position relative to c. g. (ft). 


DL.DB1 Direct flapping stiffness (rad/sec 2 ) 


DL.DA1 Off-axis flapping stiffness (rad/sec 2 ). 

D.TR Tail rotor horizontal position (ft). 

EPS.HT Downwash field on horizontal tail relative to induced 
velocity. 

F(i) Force and moment vector. 

F(l) Total x- force component (lb). 

F( 2 ) Total y-force component (lb). 

F(3) Total z-force component (Id). 

F(4) Total rolling moment component (ft-lb). 

F(5) Total pitching moment component (ft-lb). 

F( 6 ) Total yawing moment component (ft-lb). 

F( 7 ) Pitch axis flapping angle, a 1 , (rad). 

F( 8 ) Roll axis flapping angle, bj, (rad). 

FR.MR Effective frontal area of main rotor (ft 2 ). 

FR.TR Effective frontal area of tail rotor (ft 2 ). 

GAM . OM . 16 One-sixteenth the product of Lock Number and rotor 
angular rate (rad/sec). 

GV(7) Longitudinal tip-path-plane angular rate (rad/sec). 

GV(8) Lateral tip-path-plane angular rate (rad/sec). 

H.FUS Fuselage vertical position of aerodynamic center 

relative to c. g. (ft). 

H.HUB Hub vertical position relative to c. g. (ft). 

H.TR Tail rotor vertical position relative to c. g. (ft). 


GV( 7 ) 
GV( 8 ) 
H.FUS 


H.HUB 


B-2 





HP. LOSS Net power loss due to transmission, accessories, etc. 1 

(hp) • 

IS Main rotor shaft incidence (rad). 

ITB Inverse tip-path-plane lag (rad/sec). 

ITB2.0M ITB squared over OMEGA. MR (rad/sec). 

KC Flapping coupling factor. 

KIND Induced velocity factor. 

L.FUS Fuselage aerodynamic rolling moment (ft-lb). 

L.MR Main rotor rolling moment (ft-lb). 

L.TR Tail rotor rolling moment (ft-lb). 

L. VT Vertical tail rolling moment (ft-lb). 

M Vehicle mass (slug). 

M. FUS Fuselage aerodynamic pitching moment (ft-lb). 

M.HT Horizontal tail pitching moment (ft-lb). 

M.MR Main rotor pitching moment (ft-lb). 

M. WN Wing pitching moment (ft-lb). 

N. FUS Fuselage aerodynamic yawing moment (ft-lb). 

N.MR Main rotor yawing moment or torque (ft-lb). 

N.TR Tail rotor yawing moment (ft-lb). 

N.VT Vertical tail yawing moment (ft-lb). 

OMEGA. MR Main rotor angular velocity (rad/sec). 

OMEGA.TR Tail rotor angular velocity (rad/sec). 

P . CLIMB Power loss due to change in potential energy (ft- i 

lb/sec). 

P. INDUCED. MR Power loss due to main rotor induced velocity (ft- 

lb/sec). * 

i 

P. INDUCED.TR Power loss due to tail rotor induced velocity (ft- 
lb/sec ) . 

P. PARASITE Power loss due to fuselage parasite drag (ft- 

lb/sec). 

P. PROFILE. MR Power loss due to main rotor profile drag (ft- 

lb/sec). 

POWER. FUS Power loss from fuselage aerodynamic drag (ft- 

lb/sec) . 

POWER. MR Power loss from main rotor and fuselage (ft- 

lb/sec ) . 
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POWER.TR Power loss from tall rotor (ft-lb/sec). 

POWER. WN Power loss from wing Induced drag (ft-lb/sec). 

POWER . ROTOR . MR Power loss from main rotor (ft-lb/sec). 
PRESS. RATIO Pressure ratio. 

R2 One half RHO. 

RHO Air density (slug/ft 2 ). 

R.MR Main rotor radius (ft). 

R.TR Tail rotor radius (ft). 

54 Sin of roll Euler angle. 

55 Sin of pitch Euler angle. 

56 Sin of heading angle. 

ST Numerical integration step size (sec). 

TEMP. RATIO Temperature ratio. 

TWST.MR Main rotor twist (rad). 

TWST.TR Tail rotor twist (rad). 

THRUST. MR Main rotor thrust (lb). 

THRUST.TR Tail rotor thrust (lb). 

TIME Present time (sec). 

TORQUE. MR Main rotor torque (ft-lb). 

VA(1) X-axis velocity relative to airmass (ft/sec). 

VA(2) Y-axis velocity relative to airmass (ft/sec). 

VA(3) Z-axis velocity relative to airmass (ft/sec). 


VA( 4 ) 

Roll-axis 
( rad/sec ) . 

angular 

velocity 

relative 

to 

airmass 

; 

t 

VA( 5 ) 

Pitch-axis 
( rad/sec ) . 

angular 

velocity 

relative 

to 

airmass 

f 

* 

f 

VA( 6 ) 

Yaw-axis 

angular 

velocity 

relative 

to 

airmass 



( rad/sec ) . 

VB(1) X-axis inertial velocity (ft/sec). 

VB(2) Y-axis inertial velocity (ft/sec). 

VB(3) Z-axis inertial velocity (ft/sec). 

VB(4) Roll-axis inertial angular velocity (rad/sec). 
VB(5) Pitch-axis inertial angular velocity (rad/sec). 

VB(6) Yaw-axis inertial angular velocity (rad/sed). 

VE(1) X-axis velocity relative to earth (ft/sec). 






) 

t 




VE( 2 ) 

VE( 3 ) 
VE(4 ) 

VE( 5 ) 

VE( 6 ) 

VG( 1 ) 

VG( 2 ) 

VG( 3 ) 

VG( 4 ) 

VG( 5 ) 

VG( 6 ) 

VHAT.2 
VI. MR 
VI. MR. 2 
VI. TR 
VI .TR.2 
VP(i) 
VR.TR 

VB.TR 

VT 

VT . IN . FPS 

VTA 

V . TIP 

WA.FUS 

WB 

WR 

WT 

XE( 1 ) 

XE( 2 ) 

XE( 3 ) 
XE(4) 

XE( 5 ) 

XE( 6 ) 


Y-axis velocity relative to earth (ft/sec). 

Z-axis velocity relative to earth (ft/sec). 

Euler angle rate (rad/sec). 

Pitch-axis Euler angle rate (rad/sec). 

Yaw-axis Euler angle rate (rad/sec). 

X-gust component (ft/sec). 

Y-gust component (ft/sec). 

Z-gust component (ft/sec). 

Inertial roll gust (rad/sec). 

Inertial pitch-gust (rad/sec). 

Inertial yaw-gust (rad/sec). 

Intermediate variable in thrust calculations (ft^sec 2 ). 
Main rotor induced velocity (ft/sec). 

VI. MR squared. 

Tail rotor induced velocity (ft/sec). 

VI. TR squared. 

Past value of VE(i). 

Net vertical velocity relative to tail rotor blade 
(ft/sec) . 

Net vertical velocity through tail rotor actuator disk 
( ft/sec) . 

Total airspeed (kt). 

Total airspeed (ft/sec). 

Total airspeed (ft/sec). 

Main rotor tip speed (ft/sec). 

Apparent vertical velocity on fuselage (ft/sec). 

Net vertical velocity relative to rotor blade (ft/sec). 
Net vertical velocity through actuator disk (ft/sec). 
Gross weight (lb). 

X-axis position (ft). 

Y-axis position (ft). 

Z-axis position (ft). 

Roll Euler '.ngle (rad). 

Pitch Euler angle (rad). 

Heading Euler angle (rad). 
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X.FUS Fuselage x- force (lb). 

X.GRAV Gravity x-force (lb). 

X.MR Main rotor x-force (lb). 

X.HT Horizontal tail x-force (lb). 

X. WN Wing x-force (lb). 

2 

XUU.FUS Fuselage parasite drag force (ft ). 

Y. FUS Fuselage y-force (lb). 

Y.GRAV Gravity y-force (lb). 

YMIN.VT Vertical tail stall factor (ft 2 ). 

Y. MR Main rotor y-force (lb). 

Y.TRT Tail rotor y-force (lb). 

YUU. VT Vertical tail profile drag factor (ft 2 ). 

YUV. VT Vertical tail circulation lift factor (ft 2 ). 

Y. VT Vertical tail y-force (lb). 

2 

YVV.FUS Fuselage sideward drag factor (ft ). 

ZMIN.HT Horizontal tail stall factor (ft 2 ). 

ZMIN.WN Wing stall factor (ft 2 ). 

Z. FUS Fuselage z-force (lb). 

Z.GRAV Gravity z-force (lb). 

Z.HT Horizontal tail z-force (lb). 

ZUU.HT Horizontal tail profile drag factor (ft 2 ). 

ZUU.WN Wing profile drag factor (ft 2 ). 

ZUW.HT Horizontal tail circulation lift factor (ft 2 ). 

ZUW.WN Wing circulation lift factor (ft 2 ). 

2 

ZWW.FUS Fuselage quadratic drag coefficient along z-axis (ft ). 
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APPENDIX C 

DEFINITION OF PROGRAM INPOT FILE 



i 

*. 

The following describes the input format needed to define 
the math model for a specific helicopter. The specific values 
given correspond to the AH- IS example. Individual entries are 
discussed in detail in Appendix D. 


♦ ♦♦♦t*********#*******************^^^^^^^^ 

* 

* DATA FILE FOR THE AH-S HELICOPTER 

* 

*********************** ************************ 


********************** 

* 

PARAMETERS * 

» 

********************** 


(CONFIGURATION, AIRCRAFT NAME, FS.CG, WL.CG, 
1C '2, "AH- IS", 196 , 75 


WT, IX, IY, IZ, 
9000, 2593,14320,12330, 


IXZ) 

0 


(FS.HUB, WL.HUR, IS, E . MR , I.B, R. MR , A. MR, RPM.MR, CD0, 
200 • 153 * 0 i 0 , 1382, 22, 6 , 324, 0.010 

(FS.FIJS, WL.FUS, XUU.FUS, YVV.FUS, ZWW.FUS) 

200 > 65, -30 , -275, -41 


(FS.WN, 

WL.WN, 

ZUU.WN, 

ZUW.WN, ZMAX.WN, B.UN) 

200, 

65 , 

-39 , 

-161, -65, 10.75 

(FS.HT, 

WL.HT, 

ZUU.HT, 

ZUW.HT, ZMAa.HT) 

400, 

65, 

o , 

-80, -32 

(FS. VT, 

ML. VT, 

YUU. VT, 

YUV.VT, YMAX.VT) 

490, 

80, 

o , 

-62, -50 

'FS. TR, 

WL.TR, 

R.TR, A, 

■TR, SOL.TR, RPM.TR, TWST. TR) 

521.5, 

119, 

4.25, 

6, .105, 1660, 0 


B . MR , 
? *- 1 


C. MR , TWST . MR , 
2.25, -.175, 


K 1 ) 
0 


i 

f 

4 

f 

) 


1 
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APPENDIX D 

MATH MODEL MATCHING PROCESS FOR AUGUSTA A109 II HELICOPTER 


This appendix describes a minimum-complexity math model 
! ' version of the Augusta A109 II helicopter based on available 

i flight data and flight manual information. The data were 

» furnished by the Army Aerof lightdynamics Directorate in order to 

I , provide an illustration of the parameter matching procedure ana 

verification of the resulting math model. 

It is believed that the results of this modeling process 
are sufficiently good to be used as the basis of a lateral 
control handling qualities experiment such as that performed 
under this contract. 

The general procedure followed was first to quantify the 
basic math model form using engineering data pz'ovided by the 
helicopter manufacturer. The second step was to adjust 
parameters in order to match trim data from flight and to’ add 
certain nonlinear characteristics such as downwash on tail and 
fuselage. The final step was to match dynamic response cases 

adjusting rotor model parameters. Details of the matching 
procedure are presented below. 


1. Initial Quantification of Model Parameters 

The first step was to set up the main data file for the 
math model using all available engineering data. In this case a 
fairly complete array of these data were supplied by the 
manufacturer . The following paragraphs present the initial 
quantification of parameters for each of the model components and 
a short discussion of the basis for quantification. 


Loading Parameters 


FS.CG, WL.CG, WT 

132.7 in , 38.5 in , 5401 lb 


FS.CG is the location of the center of gravity in the fuselage 
reference system in inches from the zero fuselage station For 
the A109 this can be found in the flight manual but must be 
converted from millimeters. 

WL . CG" is the vertical center of gravity location in incites 
above the zero waterline. Without a specific value, this car be 
estimated as approximately at the level of the engine. However 
it is important to determine this quantity as accurately as 
possible since a significant portion of flapping stiffness (thus 
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pitch and roll damping) results from the vertical offset of the 
rotor hub from the vertical center of gravity. 

WT" is the gross weight of the aircraft in pounds. A 
representative value can be picked from the flight manual loading 
envelope diagram, however a fairly accurate weight should be 
available for any given set of flight data. 


1590 slug-ft 2 , 6761 slug-ft 2 , 6407 slug-ft 2 , 598 slug-ft 2 

IX, IY," "IZ," and "IXZ" are the moments of inertia about the 
center of gravity in the "body" or fuselage reference line axis 
system. The first three are essential to the math model. IXZ 
can be neglected but an effect can be seen in yaw respose due to 
roll axis inputs. It should be recognized that moments of 
inertia often cannot be measured accurately and can therefore be 
subject to modification in order to match flight data. For the 
A109 the value of IX was reduced from that shown above in order 
to match the primary roll damping mode. 


Main Rotor 


FS . HUB , WL.HUB, IS. E.MR, I.B, 

132.4 in, 98.2 in, .11 rad, 0.5 ft, 212 slug-ft 2 


"FS.HUB" is the fuselage reference system location of the main 
rotor hub measured aft of the zero fuselage station. 

"WL.HUB" is the corresponding waterline location of the main 
rotor hub. 

IS" is the main rotor shaft tilt forward of vertical in the 
fuselage reference system and measured in radians. 

E.MR is the geometric main rotor flapping hinge offset for an 
articulated rotor or the effective hinge offset for a rigid 
rotor. Any empirical adjustment of this parameter should be done 
with care. In general, for teetering and articulated hubs, 
variation of the roll moment of inertia is probably easier to 
justify than the geometric flapping hinge offset. 


I.B" is the flapping inertia 
flapping hinge. 


single blade about 
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R.MR, 


A. MR, RPM.MR, C D0, 


B.MR, 


18 ft, 6/rad , 385 rpm, 0.010, 4 blades 

"R.MR" is the actual main rotor diameter in ft. 

units of non-dimensional V lift lf coeff ic~ 3 1 f Pe ° f the main rotor ir ‘ 

angle of attack. Vaiues 

units^t revolutiona n per a mlnite ^Thi 1 ^ ° f th * ” ain rot in 

a flight manual. nunute. This can ordinarily be found in 

cross sectio^ vflue^of 0*010 0*012 f ° r the main rotor blad 

- C f? be adjusted in order to^f?? 6 Commonl y used, bu 
especially in hovering flight. flt power required data, 

"B.MR" 1, the number of blades in the main rotor array. 


C. MR, 


TWST.MR, K1 


1.10 ft, -.105 rad, .096 

C.MR" i s the blade chord in ft. 

••TWST.MR- is the effective blade twist in radians. 

based * centreing hLge^eomitry* 1lUh tlV h P j; toh " fla P coupling 
given for the A109, the Although the above value was 

matching pitch and rolfcross-cMpUng eHSt' 1 ' "limited "i“ 
£u.s^Ia&^_ P arameters 


FS.FUS, 
132 in. 


WL.FUS, 
38 in. 


X00.F0S YW.FUS. ZWW.FUS 

10.8 ft , -167 ft 2 _ 85 2 


FS FUG 1 

center °of IxiT 5 ^ 1 "? to the effective 

e to the main rotorhub C pLi^on. 13 WaS ^^ally 

''WL.FUS' 1 is the waterline for 

longtudinal axis. It can be adjusted^ 61 " ° f pressur e in the 
well as that of the fuselage itsilf ^°K? CC ° Unt for hub drag as 
the vertical center of gravity. ' 1 thlS case ib was set at. 
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dra^' F in ef f ectiv ® f ront f 1 area corresponding to profile 

manufacturer x " axis - The value used here was provided by the 

'YVV.FUS" is the effective side area for sideward fliaht that 

tD 90 d9g - A8al " “ - providid t; 

5ur^ e fiJspe2nh?nges 3 U8ed in conjuntiori with Pitching moment 


FS.WN, WL.WN, ZUU.WN, ZUW.WN, ZMAX.WN, B.WN 
000 in, 00 in, 000 ft 2 , 0000 ft 2 , 000 ft 2 , 1 ft 

The A109 does not have a wing, thus zeros were set for all valu** 
except the span which needs any arbitrary non-zero value to avoid 

^ i!i° n ^ 7 H Z !y°: , the indiv idual values can however be found 

or estimated similarly to those for the horizontal tail. 




FS.HT, 


WL.HT, ZUU.HT, ZUW.HT, 


ZMAX.HT 


330 in, 54 in, .4 ft 2 , -34 ft 2 , -22 ft J 


tail^in^ ttve aerodynamic center of the horizontal 

estimated fr ° m ^ reference fuselage station. It can be 

enairlivf „ ? + the quarter mean aerodynamic chord based on 
engineering data or on a planview of the aircraft. 

t.aii HT TL th v,? f ,! eCti a e - vertical location of the horizontal 
• ii. The value used is important in computing the position of 
the mam rotor downwash field as airspeed is va?fed PO:3ltlon of 

ana?» H If t Jj® e ? f J ctive Hit per unit dynamic pressure at zero 

! °f “ tack P®l®tlve to the fuselage referenci system. The 

high^forward vIloc^tTe^' ‘ B ® 3tabllshi >'* trim pitch angle at 

ZUW.HT is the effective variation in circulation lift and can 
surface m arJa. aS th ® negative Product of lift-curve-slope and 


D-4 


RP*vr 




FS.VT, WL.VT, YUU.VT, YUV.VT, YMAX.VT 

380 in, 80 in, 3.3 ft 2 , -47 ft 2 , -17 ft 2 

"FS.VT'' is the fuselage station for the effective aerodynamic 
center of the vertical fin. 

"WL.VT" is the vertical position of the vertical fin aerodynamic 
center . 

"YUU.VT" is the net y-force per unit dynamic pressure for zero 
sideslip. This arises either from vertical fin camber or 
incidence . 

"YUV.VT" is the sideforce arising from a side-velocity component 
and is approximately equal to the lift-curve slope time the net 
fin area. 

YMAX.VT sets the maximum sideforce generated by the vertical 
tail at stall. 


FS.TR, WL.TR, R.TR, A.TR, SOL.TR, RFM.TR, TWST.TR 

391 in, 70 in, 3.1 ft, 3/rad, .134, 2080 rpm, -.137 rad 


"FS.TR" and "WL.TR" represent the center of the tail rotor hub in 
the fuselage reference system. 

"R.TR" is the radius of the tail rotor in ft. 

"A.TR" is the effective lift-curve-slope of the tail rotor and 
can be set equal to that of the main rotor. It can be adjusted 
downward in order to account for interference effects with the 
vertical fin. In this case it was reduced by one half in order 
to match pedal trim data as discussed below. 

"SOL.TR" is the solidity of the tail rotor, i. e. , the ratio of 
actual blade area to disk area. 

"RPM. TP." is the angular velocity of the tail rotor in terms of 
revolutions per minute. 

"TWST.TR" is the effective twist of the tail rotor blade. 


/ 




2. Model Adjustments Needed to Mutch Static Trim Data 

Several model parameters and functions were adjusted in 
order to produce good static trim matches. This was important in 
establishing realistic attitudes, control deflections, and power 
requirements. 

For the A109 those adjustments found necessary included: 

° A tail rotor efficiency factor of 0.5 in order to account 
for vertical fin interference effects. 

o A triangular main rotor downwash field superimposed on the 
horizontal tail vertical velocity and displaced 1 ft 
rearward. 

o A magnification of the dihedral effect at low speeds 
individually set for the lateral and longitudinal axes. 

o A shift in the planview center of pressure such that the 
downwash on the fuselage provides a pitching moment 
proportional to airspeed. 


Tail 


The first adjustment was made by changing "A.TR" from a 
nominal value of 6/rad to 3/rad. This was done on the basis of 

matching the pedal deflection (i. e. , tail rotor collective 

pitch), especially at low speed and hover. The effect was 

applied directly in the tail rotor thrust equation wherein "A TR” 

appears. 

THRUST. TR .= (VB.TR - VI. TR) 

♦OMEGA. TR*R. TR^RHO^A. TR^SOL. TR^PI*R. TR^R. TR/4 


. ^he secon d adjustment consisted of assuming a triangular 
induced velocity field with a magnitude of 2 at the rotor tip and 
aero at the hub. This can be justified by observing measured 

f i e il^ ata such 35 Presented by Heyson and Katzoff in 
NACA Report 1319. This effect is crucial to portraying the large 
c ange in pitch attitude between zero airspeed and 20 kt 
rearward. The program instructions affected are given below: 

D.DW=( VA ( 1 )/( VI .MR-VA^ 3 )*(H.HUB-H.HT) )-( D. HT-D. HUB-R.MR ) 

horizontal tail)* 6 6dge ° f downwaah P asses through plane of 
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D.DW=D.DW + 1 

(shift of D.DW by one foot in order to match speed where downwash 
on tail effect is seen in trim data) 

IF (D.DW>0 AND D. DW<R.MR) THEN EPS.HT= 2*( 1-D. DW/R.MR) 

ELSE EPS . HT-0 

(triangular downwash if D.DW is negative) 

WA.HT ~ VA(3) - EPS. HT*VI .MR .+ D.HT*VA(5) :* local z-vel at h.t. 
(appearance of downwash effect in computation of relative z- 
velocity component at horizontal tail) 


Dihedral Magnification at Low Sneed 

The third adjustment consisted of magnifying the effective 
dihedral effect at very low speeds when the rotor wake interacts 
with the fuselage. This enhanced dihedral effect could be seen 
directly in the cyclic control gradient with respect to forward 
speed and side velocity. 

The model equations affected are limited to the rotor 
flapping equations. Below a speed of VTRANS the computed 
parameters daldu and dbldv are multiplied by 3 and 2, 
respectively. The values are empirical and based on cyclic trim 
data. VTRANS was set at 30 kt based on the large change in stick 
trim observed at that point. The program statments involved are 
shown below: 

IF VA( 1) < VTRANS THEN WAKE.FN = 1 ELSE WAKE. FN = 0 

(rotor wake effects are added to the effective tip-path-plane 
dihedral when WAKE.FN = 1, i. e., below and airspeed equal to 

VTRANS) 

A. SOM = GV(8 )-DC( 2)+KC*GV(7 )+DBlDV*VA( 2 )*( 1+WAKE. FN) 

( i. e., bl - Al + e.al + dbl/dv .V ) 

B. SUM = GV(7)+DC(3)-KC*GV(8)+DA1DU*VA(1)*(1+2*WAKE.FN) 

( i. e. , al + Bl - e.bl + dal/du .0 ) ; 

GR(7)= - ITB*B. SOM ITB2.0M*A. SOM - VA(5) 

( i . e . , al . dot = . . . ) 

GR( 8 ) = - ITB*A. SUM + ITB2 . OM*B . SOM - VA(4) 

( i . e . , bl . dot = . . . ) 


Ed-vmwash Center of Pressure on Fuselage 

The fourth and final adjustment needed to match trim data 
is the shift of downwash center of pressure on the fuselage as 
speed varies. This affects not only the longitudinal cyclic to 
trim but also the trim pitch attitude. 




The approach was to compute an effective wake postion in 
the plane of the fuselage similar to that computed for the 
horizontal tail. This position was then used in the pitching 
moment equation along with an empirical magnification factor. 
The net effect is a change in fuselage pitching moment with 
forward speed. 

The program instructions affected are: 


WA.FUS = VA( 3 ) - VI. MR 

(computed net downwash on fuselage, i. e. , W - V ) 

D. FW= ( VA( 1 )/( -WA. FUS )*(H. HUB-H. FUS) )-( D.FUS-D.HUB ) 

(computed position of downwash at fuselage waterline as airspeed 
varies ) 

D.FW=3*D.FW 

(empirical magnification of a.c. shift used to match trim data) 

Z.FUS - R2 * ZWW.FUS * ABS(WA.FUS) * WA.FUS 

(s-force resulting from downwash on fuselage) 

M. FUS = Z.FUS * D.FW - X.FUS * H.FUS 

(pitching moment due to x- and s-forces acting at their 
respective aerodynamic centers) 

The resulting math model trim characteristics are 
compared with the A109 flight data in the following pages. 
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3. Model Adjustments Needed to Match Dynamic Response Data. 


, . ^ ie f f na l step in the model matching process was 
adjustment in the model to account for features seen in both 
primary on-axis response as well as cross-coupling effects. 


the 

the 


The adjustment of dynamic response 
limited only to pitch and roll response in hover, 
flight conditions would be addressed in a similar 


features was 
Other axes and 
manner . 


In the case of the A109 it was found that 
cross-coupling in flapping and hub moment equations 
good match could be achieved for the pitch and roll 
only minor adjustment of moments of intertia 


by removing 
a reasonably 
axes with 


fell. In ertia Reduction 

„ dominant roll-damping mode was matched closely bv 

varying the lateral flapping stiffness via a reduction in the 
roll moment o f Intertia. This was considered preferable to 

hinge offset since the latter would also 
affect the pitch response. 


The value of IX in the data input file was reduced 
1590 to 1300 slug-ft 2 . 


f ruin 


fia.iHS.yal of Rotor Flap Crosa-Omnl j n F 

attributed ° f pitGl \ res P ons e due to a roll input can be 
In t K' f + croas -coup 1 i ng in the rotor flapping equations. 
In this case it was found that decoupled flapping provided - 
better match to flight data thus the simplification was made 


The coupled first order flap equations 
simply by recomputing the values of " ITB'' and 
the primary flapping response consists only of 
as described in Reference 1 . 


were decoup Led by 
"ITB2.0M.'" Thus 
a first-order lag 


This change is accomplished by setting: 

ITB = GAM.OM.16 


and ITB2.0M = 0 


Eli.nnnat i sn — sf — Hub Moment Q ross-Coupi ing 
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The aerodynamic cross-coupling represented by DL.DA1 was 
also set equal to zero in order to further suppress pitch cross - 
coupling due to roll as seen in flight data. 


Removal of Delta- 3 Effect 

In order to maintain consistency of model complexity 
following the above simplifications, the Delta-3 effect as 
represented by the parameter "Kl" was also set to zero. 


Adjustment of Cross-Axis Inertia 

The effect of an inclined principal axis of inertia could 
be seen readily in the short-term yaw response following a roll 

input. An increase in I from 598 to 800 slug-ft 2 provided a 

X z 

slightly better match to flight data. 


Adj ustment of Cyclic Control Phasing 

There was some evidence of cyclic control phasing in the 
A109 although control system geometry indicated none. First, 
inspection of swashplate angle records for pitch and roll inputs 
showed minor off-axis inputs which were generally consistent with 
the cross-coupling response which followed. Also there was a 
direct measurement of a nearly 10 deg steady lead-lag component 
which, depending upon the hub pitch control geometry, could 
contribute to a cyclic control phase effect. 


The results of the dynamic response adjustments for pitch 
and roll inputs in hover are shown in the following pages. 


i 
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4 Revised Program Listings for the A 109 Math Model 

... The following listings show the revisions made to the data 
rile, preliminary calculations, and the dynamics subroutine. 


type HEL302.DAT 

# 

* DATA FILE FOR THE A109 II HELICOPTER PARAMETERS # 

* » 
HHtiHIHIHMHHtlHHHHHNHHHHHHIHHIHHIiHHIfHtH 

(CONFIGURATION, AIRCRAFT NAHE, FS.C6, NL.C6, NT, IX, IV, 12, HZ) 

302, -A109 IP, 132.7 , 38.5 , 5401, 1300, 6760, 6407, 800 

(FS. HUB, UL. HUB, IS, E.MR, I. 8, R.HR, A. HR, RPH.HR, CD0, 8. HR, C.HR.TNST.HR, Kl) 
132.4, 98.2, .11, 0.50, 212, 18, 6 , 385, 0.010, 4, 1.10, -.105, 0 

IFS.FUS, NL.FUS, XUU.FUS, YVV.FUS, ZHW.FUS) 

132, 38, -10.8, -167, -85 

(FS. HN, NL.NN, ZUU.NN, ZUN.NN, ZHAl.NN, B.HN) 

0 1 0 , 0 , 0 , 0 , 1 

(FS. HT , NL.HT, ZUU.HT, ZUN.HT, ZHAX.HT) 

330, 54, .4, -34, -22 

(FS. VT , NL.VT, VUU.VT, YUV.VT, YHAI.VT) 

380, 80, 3.3, -47, -17 

(FS. TR, NL.TR, R.TR, A.TR, S0L.TR, RPH.TR, TNST.TR) 

3 ’». 70, 3.1, 3, .134, 2080, -.137 
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1 490 ' If Iff If ff 

1500 ' 

1510 ' Preliminary Calculations 

1520 ' 

1S30 ' f If f fllfflfflff f fff ftff ffffffff If If tfftlf ff fill If Hff If f f fill Iff f If f Iff f 

1540 1 

1550 VT.IN.FPS * ABS ( XDOT /COS ! 6AHIM.RA0) I 
1540 VT * VT. IN. FPS/FPS. PER. KNOT 

1570 VT.IN.FPS. SQUARED * VT.IN.FPS A 2 

1560 H * NT/6RAV ! 

1590 ONESA.HR = RPH.HR«2fPl/40 

1600 OHESA.TR = RPH.TRf2fPI/60 

1605 V.TIP*R.HRfQHE6A.HR 

1610 FR.HR = CD0*R. HRfB.HRfC.nR 

1620 FR.TR * CDOfR. TRfB.TRfC.TR 

1630 . HP. LOSS = 90 , 

1640 VTRANS * 50 speed lor transition Iron dihedral wake function I 

1660 TEHP. RAT 10=1 ! - LAPSE. RTtH 

1670 PRESS. RATIO = TEHP.RAT10 A TEHP.EX? 

1660 DENS. RATIO = PRESS. RATIO/TEHP. RATIO 

1665 RO = DENS. RAT IOfRHO. SEA. LEVEL : R2*R0/2 

1666 6AH.0H.16 = ROf A. HRfC. HRtR. HR A 4/ I . B fOHEGA.HR/16f < l+8/3f E. HR/R.HR) 

1690 KC = ( .75fGNE6A.HRfE.HR/R. HR / BAH. OH. 16 ) + Ki: ' flapping aero cpl 

1691 1TB2.0H * 0HE6A. HR/ (1* I0HE6A.HR/GAH. OH. 16) A 2) s' flapping x-cpl coef 

1692 1TB - 1TB2. OHf CHE6 A.hr/GAH. OH. 16 flapping prinary resp 

1693 FTB2.0H = 0 : ITS = GAH.0H.16 teap nod for A109 * 

1694 DL.DB1 = B. HR/2< 1 1 . 5f I . BfE . HR/R. HRfOHEGA.HR* OHEGA. HR) prinary flapping stiffness j 


1696 

0L.DAUR2fA.HRfB.HRfC.HRfR.HRfV.TlPfV.TIPfE.HR/6:' cross flapping stiffness 

1697 

CT=WT/ (ROf PUR.HRf R. HRfV. T1P«V.TIP) 

s’ thrust coefficient 

1696 

A. SI 6HA=A. HRf B. HR«C. HR/R-HR/PI 

: ' a x signa 

1699 

DB1DV=2/0HEGA. HR/R. HRf (8»CT/A. SIBHA+ISQRICT/2) ) ) : ' TPP dihedral effect 

1700 

DA1DU=-DB1DV 


: ' TPP pitchup with speed 

1702 1 





1730 

H.HUB = (NL.HUB-HL.CGI/12 

D.HUB * 

IFS.HUB-FS.C6I/12 

' hub re eg 

1740 

H.FUS = IHL.FUS-NL.C61/12 

D.FUS - 

(FS.FUS-FS.C6I/12 

' fuselage re eg 

1750 

H.HN - (HL.HN -HL.CG1/12 

3.HN * 

(FS.HN -FS. CB) / 12 

' wing re eg 

1760 

H.HT * (NL.HT -NL.C61/12 

D.HT * 

(FS.HT -FS. CGI/12 

' horizontal tail re eg 

1770 

H.VT * IHL.VT -NL.C61/12 

D.VT = 

IFS.VT -FS.CG1/12 

' vertical fin re eg 

1780 

H.TR * 1NL.TR -HL.CG1/12 

O.TR * 

(FS.TR -FS.C61/12 

' tail rotor re eg 

1790 ’ 





1800 

RETURN 
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Dynamics Subroutine 


ORIGINAL PAGE IS 
OF POOR QUALITY 


6000 *#mmi HHHHHHHHNHHHNHHHIHtHHHHHHHIiHHHHHi 
0010 ■ 

6020 ' DYNAMICS: Dysaaics subroutine 

6030 ‘ 

6040 ' ttlHIiHHIiHHHHHHHHIHHHIHHHHHIHiHIIHHIiHMHHtH 
6070 ‘ 


6080 
6090 ' 
6100 
6120 ' 
6130 
6140 
6160 1 
6170 
6180 
6190 
6200 
6210 
6220 
6230 ‘ 
6290 
6320 * 
6330 ‘ 

6340 ‘ 

6341 

6345 

6346 

6347 ’ 
4350 
6360 
6370 ' 
5380 
6390 
6400 ‘ 
6410 
6420 ' 
6430 ' 
6440 ' 

i . 6490 

fc 6500 

6510 ' 
6540 
6550 
65/0 
6580 
6590 
6610 
6620 * 


IF CHECK * 1 THEN 60SUB 11940 : ' print state variable to screen 

(HHHHHHit Prelieinary calculations mhhhhh« 

C4 * C0SUE(4I) : S4 * SINCXEC4)) evaluate Euler angle trig 4ns 
C5 * COSUE (5)) : SS * SINIXEC5D 

VA(l)sVBIl) -<V6(l)eC5> evaluate relative airaass velocities 
VA(2)*VB(2) *1V8(2I*C6-V6(1)«S6) 

VAI3)=V8(3) :* -IV6(3)*C5*»GI1)«S5) 

VAI4)*VB(4) -V614) 

VA(5)«VB(5» 

VA(6)*VBI6) 


VTA*SBR(VAIl)eVA(l)*VA(2)eVA(2)+VA(3)eVA(3)) 


hhhhhhhh Rotor tip path plane dynaaics eteeeeeeeeeeee 


IF VAC 1 1 < VTRANS THEN NAKE.FN * 1 ELSE NAKE.FN = 0 : ‘ rotor uake eCfects 

A.SUH« BV( 8 )-DCI 2 l*KCf 6 » 17 )*DBlDV*VA« 2 )«(l+NAKE.FN) bl - A1 ♦ e.al ♦ dbl/dv .V 

8 .SUH« 6 V ( 7 ) +DC ( 3 ) -KC* 6 V ( 8 ) +DA 1 IHHVA ( 1 ) # 1 1 + 2 sHAKE. FN) al ♦ Bl - e.bl ♦ dal/du .U 


6R<7>* - ITBeB.SUH - ITB2.0HeA.SUH - VACS) al.dot 

6R(81* - ITBsA.SUH ♦ ITB2.0HeB.SUH - VAC4) bl.dot 


8VC7)*6VC7) ♦ STe(A2»6R(7) ♦ B2«APC7U al updated 

6V(B)«6V(8> ♦ STe(A2efiR(8) ♦ B2eAP<8H bl updated 


AF(7)«6RI7> : APC8)*6RC8) 


save past values 


eteeeeeeeeeeee Main Rotor thrust and induced velocity eeeeeeeeeeeee 

HR * VAC3) ♦ C6V C7> - IS)eVA(l) - 6VC8)eVAC2) r z-axis vel re rotor plane 
N8 • NR +2/3eOHE6A.HReR.HRe(OC(l) ♦ .75eTNST.HR)i':-aitis vel re blade 

FOR 1*1 TO 5 iterative solution o( thrust and induced vel 

THRUST. HR*(NB-VI.HR)eOHE6A.HReR.HR«RHOeA.HR«B.HRtC.HReR.HR/4 
VHAT . 2*VA C 1 ) *2 ♦ VAI2) A 2 ♦ NRe(NR*2«Vl.HR) 

VI . HR. 2*S8R I (VHAT. 2/2)e (VHAT . 2/2) ♦ (THRUST. HR/ 2/ (RHOePle R. NR A 2) ) A 2) - VHAT.2/2 
VI.HR*S0RIABS(VI.HR.2)I asm rotor induced velocity 
NEXT 1 


r 




t 


A 
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Fuseli?* HIHHHIHHHtHHtllMHH 


6630 • 

6640 ' 

6650 NA.FUS * VA(3) - VI. HR include rotor dounuash on fuselage 

6652 D.FN*( VA( 1 ) / (-HA.FUS) * CH. HUB-H. FUS) M 0.FUS-D.HU9 Is' position of dowiMsh on fuselage 

6653 D.FN*3«D.FN : ' empirical correction for fus a.c. shift eagnification 
6660 ' 

6670 I. FUS « R2 • KUU.FUS * ABS(VAU) ) • VA(l) :* drag force 

6680 Y.FUS * R2 * YW.FUS * ABS(VA(2)) • WH2) tide-force 

6690 2. FUS * 82 * JNN.FUS i ABS (NA.FUS) « HA.FUS heave force 

6710 L.FUS * Y.FUStH.FUS 

6720 H. FUS * Z.FUS«D.FN - l.FUStH.FUS 

6740 * 

6780 F. 1N8UCE0.NR » THRUST. HR « VI. HR 
6790 P.CLIH8 « NT»HD0T 

6800 F. PARASITE * - I.FUS#VA(1> - Y.FUS»VA<2) - 7.FUS*M.FUS 

6810 P. PROFILE. HR * R2* (FR. HR/4) <0HE6A. HR*R. HR* (ONEGA. HR A 2<R. HR A 2 ♦ 4.6*(VA(1>«VAU)+VA(2)«VA(2») 
6820 PONER.HR * P.INOUCED.HR ♦ P.CL1HB ♦ P.PARASITE ♦ P.PROFILE.HR 
6830 PONER.ROTOR.HR * P.INDUCED.HR ♦ P.PROFILE.HR 

6840 POWER. FUS * P. PARASITE 

6850 TORBUE.HR = P0NER.HR/0HE6A.HR 

6860 ‘ 

6870 ‘ Coapute oain rotor force and eoaent coepocents. 

6860 ' 

6890 I.HR * -THRUST.HR * (GV(7)-IS) 

6900 Y.flR = THRUST.HR # 6V(B> 

6910 2. HR = -THRUST.HR 

6960 L.HR * Y.HRtH.HUB + DL.DBH6V18) * DL.DA1«(6V(7HDC(3) - KH6V<6!) 

6970 H.HR * 2.HR*D.HUB - I.HR*H.HUB ♦ DL.DB1»GV(7) ♦DL.DA1»(-6V(B(*DC(2) - K1*GV(7)) 

6980 N.HR * TORQUE.HR 

6990 ' 

7000 ' hhhhihih Tail Rotor thrust and induced velocity efttemtetf* 

7010 ' 

7060 VR.TR ~ *(VA(2) - VAI6MD.TR ♦ VA f 4) *M. TR1 velocity relative to rotor plane 

7070 VB.TR = VR.TR +2/3t0HE6A.TR*R.TRe(DC(4)+TNST.TR«.75) s' velocity relative to blade 
7100 - 

7110 FOR 1*1 TO 5 iterate on thrust and induced velocity 

7120 THRUST . TR* l VB. TR-VI . TR) »0HE6A. TR*R. TR*RH0*A. TR«S0L. TR»PI»R. TR»R. TR/4 

7U0 VHAT.2*(VA(J)+VA(5)«D.TR) A 2 ♦ VA(1) A 2 ♦ VR.TRe(VR.TR-2«VI.TRl 

7150 Vl.TR.2*SQR((VHAT.2/2)«IVHAT.2/2)»(THRUST.TR/2/(RH0«Pl»R.TR A 2)) A 2) - VHAT.2/2 

7160 VI. TR*S6R (ABSIVI.TR. 2)1 

7180 NEXT 1 

7190 ' 

7220 PONER.TR * THRUST. TReVI.TR 

729C Y. TR * THRUST.TR 

7310 L.TR * Y.TR*H.TR 

733C N.TF. * -Y.TR»D.TR 

7240 ' 


7360 

7370 

7410 

7411 

7415 

7416 
7440 
7450 
7455 
7460 
7470 
7500 
7530 


•HHiHHimimM Horizontal tail mmmemmatfememne 

D.DN*( VA tl ) / (V!.HR-VA*3)* IH. HUB-H. HT) )-( D.HT-D.HUB-R.HR Is' iepsngeeent of dovnvash on tail? 
0.DN*D.IN ♦ ! s shifts transition position for eepirical correction 
EPS.HT*.5»(!+S6N(D.DN)) s' uni fort dounuash field 

IF (B.Oh'O AND D.DN<R.HR> THEN EPS.HT* 2»(I-D.8N/R.HR) ELSE EPS.HT-0 s' triangular dounuash field 

NA.HT * VAI3) • tPS.HTeVI.HR ♦ D.HT»VA<5) s' local i-vel at h.t. 
VTA.H(«SQfi(VA(lltVA(ll*VA(2)«VA(2)*NA.HTeNA.HT) 

Z.HT»A2«(2UU.HT«ABSIVA(I>)MI1) ♦ 2UN.HT»ABS(VA(1))»NA.HT) s' circulation lift on h.t. 


IF ABS(NA.HT)>,3*ABS(VA(1I) THEN 2 . HT*R2«2HAI .HToABS I VTA. HT ) »NA. HT s' surface stalled’ 
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7560 H.HT * 2.NT«D.KT pitching aoeent 
757: 

7590 ' HHHHHHHHHimHH Ming HtWHIlIttHIHIIKimmiHIti 

7600 

7640 HA. BN * VAI3) -VI. NR : local z-vel at Ming 
7650 VTA. BN =SeS(VAIl>»VA(l)4BA.NNiNA.NN! 

7660 • 

7670 2 . BN*32» I ZUU. BNt VA II ) t VA ( I ) 4 2UB.NN<VAU)»HA.BN) s 1 nor.il force 

7660 1. BN--R2/PI/VTA. BN/VTA. BB4(2UU.BN«VA(1)»VA(1)42UB.BN#VA(1)»BA. BN) A 2 : induce Brag 

7690 • 

7720 IF A8S(NA.NN)>.3aABSIVA(]» THEN 2.NN*R2»ZHAI.BN4ABSIVTA.NN)4NA.NN surface stalled? 

7750 • 

7820 POBER.BN * ABS(I.BN«VA(1>> 

7840 POWER * P0BER.HR ♦ P0BER.TR 4 POBER.BN 4 HP.L0SS*55O 
7850 1 

7860 ' HHHHHHI<I<HHH< Vertical tail Mteestseeseemeestetetem* 

7670 ‘ 

7880 ' Coapute aerodynaeic forces on vertical tail 

7890 

7895 VA.VT=VA(2)4V1.TR-D.VT»VA(6) 

7897 VTA. VT*S0R<VA ( 1 ) »VA( 1 ) 4VA. VT»VA. VT) 

7900 Y. VT «R2* (YOU. VT*ABS ( VA 1 1 ) 1 »VA ( 1 ) 4 YUV.VT«ABS(VA(1))*VA.VT) 

7910 

7940 IF kBS(VA.VT)>. 3»AB5(VAI1)) THEN Y.VT*R2«YHAl.VTeABS<VTA.VT)»VA.VT surface stalled’ 

7970 ' 

8000 L.VT = Y.VTtH.VT 

8010 N.VT * -Y.VTeD.VT 

8020 

8030 mt»«»e#«a#f»ee 6eneral force equations 

8040 ' 

8050 I.6RAV - -H»BRAV«S5 gravity forces 

8060 Y.6RAV * H»GRAV»S4»C5 

8070 2.6RAV * H»6RAV«C5»C4 


8090 

8100 

8110 

8120 

8130 

8140 

8150 

8160 

8170 


8190 

8192 

8193 
8200 
6210 
8270 
6290 
8300 
8310 
8320 
833* 
8340 
8350 
8360 
8370 
8380 



gravity 

1 

H.R. 

FOS. 

T.R 

H.T. 

N1N6 


i 







F(l) * 

1.6RAV 4 I.HR 

4 l.FUS 4 

I.ttN 

F (2) * 

Y.BRAV 4 Y.HR 

4 Y.FUS 4 Y.TR 

♦ 

F (3) = 

2.6RAV 4 2. HR 

4 l.FUS 4 ♦ 2.HT 4 

2. KM 

F (41 = 

4 l.HR 

4 L.FUS * L.TR 

♦ 

FIS) * 

4 H.HR 

4 H.FUS 4 H.HT 


F(6) - 

4 N.HF: 

4 N.TR 

♦ 


V.T. ! cocponent 


Y.VT 

L.VT 

N.VT 


I-force 

Y-force 

7-force 

L-aoacnt 

H-eoeent 

N-eoaent 


F(7) 

FIB) 


BR(7)/!TB 

GR!B)/ITB 


pitch flap 
roll flap 


IF CHECKS THEN 60SUB 11040 fill force coaponent array 


Body Accelerations 

ABU) * - (VB(5)«VBI3)-V8I4HVBI2>) 4 FUJ/H 

AB(2) » (VBI4)*VBI3)-V8I1)*VBI6>) 4 FI2)/H 

AB(3) * (VBll)tVB<5)-VBI4)*VB(21) 4 FI3)/H 

A8i4) * F(4)/II 

AB(5) * FI51/1Y - VBI4)#VB(6)t(lM2)/IY 4 (VB(6)*VB(6)-VB(4)#VB(4) )«1»/IY 

AB(6) > FI61/12 4 II2*AB(4)/12 
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Integrate Body Accelerations 


FDA II * 1 TS 6 

VltlX) « VIIIX) 4 ST e (Ai e Mill) 4 |i e AP(ID) 
AP(1X) > AB(Il) : REA SAVE ACCEL PAST VALUES 
NEIT IX 


Translore to earth IA/C rel to deck) velocities 


VE(1> * (Vt(l) « C5 4 VI (3) • S3) e C4 • COS (IE(4)I 
VE(2» * V0(2)*C0S(IE(A))4VI(1) e SIR (11(A)) 

VE(3> * (VOID e S5 - VII3) • C5 ) * C4 

VEI4) * VK4) 4 (VI (S) « S4 4 VI (4) e C4) e TAR(IE(S)I 

VE IS) * VBIS) * C4 - VK6) • S4 

VE(4) • (VB16I « C4 4 VI (3) * S4> / C3 


Integrate earth (A/C relative to deck) velocities 


FOR II * 1 TO 6 

XE(1X) - IE (IX) 4 ST e (A2 » VEIIX) 4 12 e VP(IZ)) 
VP(IX) * VE(IX) : REH SAVE VEL PAST VALUES 
NEXT IX 


T!HE*T!R£4ST 


IF CHECK 1 1 THEN IF CHECK. LOOPCCHECK. LOOP. RAX THEN 60T0 4020 
IF CHECKS THEN 60SUB 12140 


RETURN 


ORIGINAL PAGE IS 
OF POOR QUALITY 


D-24 


J 


2. Gouammtnt Acomion No. 


Report Documentation Page 

S'-Cf A<*iwi*rjfa,i W 


1. Report No. 

USA AVSCOM TR-87-A-7 
NASA CR 177476 


4. Title and Subtitle 

Minimum-Complexity Helicopter Simulation Math Model 


3. Recipient'* Catalog No. 


5 Report Data 

April 1988 


6. Performing Organization Code 


7. Authoris) 

Robert K. Heffley 
Marc A. Mnich 


9. Performing Organization Name and Addreea 

Manudyne Systems, Inc, 

349 First Street 

Los Altos, California 94022 


12. Sponsoring Agency Name and Addreee 
National Aeronautics and Space Administration, 
Washington, D.C. 20546-0001, and US Army Aviation 
Systems Command, St. Louis, MO 63120-1798 


8. reforming Organization Report No 


10. Work Unit No. 

992-21-01 


11. Contract or Grant No. 

NAS2-1 1665 


13. ’vp, of Report wtd Period Covered 

Final Contractor Report 
7/AS to 7/A7 


14. Sponsoring Agency Code 


IS. Supplementary Note, 

Point of Contact: 

Michelle M. Eshow (Contract Technical Monitor) 

Aerof lightdynamlcs Directorate 

M/S 211-2, Ames Research Center, Moffett Field, CA 94035 
(415) 694-5272 FTS 464-5272 


16. Abstract 

An example of a minimal complexity simulation helicopter math model is presented. 
Motivating factors are the computational delays, cost, and inflexibility of the 
very sophisticated math models now in common use. A helicopter model form is 
given which addresses each of these factors and provides better engineering 
understanding of the specific handling qualities features which are apparent to the 
simulator pilot. The technical approach begins with specification of features which 
are to be modeled followed by a build-up of individual vehicle components and 
definition of equations. Model matching and estimation procedures are given 
which enable the modeling of specific helicopters from basic data sources such as 
flight manuals. Checkout procedures are given which provide for total model 
validation. A number of possible model extensions and refinements are discussed. 
Math model computer programs are defined and listed. 


17. Key Word* (Suggested by Authorial) 

helicopter, helicopter simulation, mathe- 
matical modeling, real-time simulation, 
handling qualities, helicopter dynamlca 


IS. Distribution Statement 


Unclassified-Unlimited 
Subject Category 08 


IS. Security Claeatf. (of thia report) 

20. Security Claaalf. (of thia page) 

21. No. of pggn 

Unclassified 

Unclassified 

100 





















